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IMPROVEMENTS TO THE MISSILE AERODYNAMICS 
PREDICTION CODE DEMON3 


Mamix F. E. Dillenius 
David L. Johnson 
Daniel J. Lesieutre 

Nielsen Engineering & Research, Inc. 


SUMMARY 


Computer program DEMON3 was developed for the aerodynamic analysis of 
nonconventional supersonic configurations comprising a body with noncircular cross 
section and up to two wing or fin sections. Within a wing or fin section, the lifting surfaces 
may be in cruciform, triform, planar, or low profile layouts, and the planforms of the lifting 
surfaces allow for breaks in sweep. The body and the fin sections are modeled by triplet 
and constant u-velocity panels, respectively, accounting for mutual body-fin interference. 
Fin thickness effects are included through the use of supersonic planar source panels. One 
of the unique features of DEMON3 is the modeling of high angle of attack vortical effects 
associated with the lifting surfaces and the body. In addition, shock expansion and 
Newtonian pressure calculation methods can be optionally engaged. These two- 
dimensional nonlinear methods are augmented by aerodynamic interference determined 
from the linear panel methods. Depending on geometric details of the body, the DEMON3 
program can be used to analyze nonconventional configurations at angles of attack up to 25 
deg for Mach numbers from 1.1 to 6. Calculative results and comparisons with 
experimental data demonstrate the capabilities of DEMON3. Limitations and deficiencies 
are listed. 


INTRODUCTION 


This final technical report describes the work performed in accordance with Tasks 1, 2, 
6, 7, and 8 specified in the Statement of Work under Section C of Contract NAS1-17077, 
Amendment/Modification No. 17. These tasks are associated directly with the engineering 
level DEMON3 computer program designed for the aerodynamic analysis of supersonic 
configurations consisting of a forward fin section and a tail fin section mounted on a body 
with arbitrary cross section. The initial development of the DEMON3 program was started 
under Contract No. 1-16469. As a result of the initial findings, certain fundamental 
problems related to the supersonic panel solutions were defined. Additional tasks to 
alleviate the problems, to add nonlinear effects of body vorticity, and to extend the range of 
applicability in terms of Mach number were included in the statement of work referenced 
above. Specifically, the Contractor performed work to accomplish the following tasks. 

Task 1. Improve the fin-on-body interference account by the incorporation of 
constant u-velocity panels in the interference shell; retain triplet panels to 



model noncircular cross section body-volume and angle-of-attack effects. 
Apply the improved program to sample cases and to cases for which 
experimental data are available. 

Task 2. Automate the sequential use of the subprograms by the incorporation of 
the driver routine designed to manage the step-wise procedure and 
exchange of data sets including fin vortex characteristics between the 
subprograms; to the maximum extent possible, simplify program input. 

Task 6. Update DEMON3 as follows: Extend fin pressure calculations to include 
as an option the combined nonlinear/ linear pressure calculation methods; 
make arrangements for exchanging data sets with NOSEVTX (see below); 
improve fin thickness routines to handle fins that are not located in radial 
planes. This task continues and modifies the effort described for Task 1 
above. 

Task 7. Join DEMON3 and NOSEVTX to treat a complete configuration consisting 
of a forward-finned section and a tail-finned section mounted on a body 
with arbitrary cross sections. Develop an executive routine to manage the 
stepwise procedure and to organize the exchange of information between 
module DEMON3 and NOSEVTX. Make comparisons with experimental 
data. This task will continue and modify the efforts described in Task 2 
above. 

Task 8. Validate the methodology developed under the Tasks 1, 2, 6, and 7 by 
comparison with experimental data. Describe the methodology and 
results in a final technical report. 

The DEMON3 computer program is an improved, extended, and unified version of the 
earlier DEMON2 set of separate programs for analyzing supersonic monoplane or cruciform 
wing-body-tail combinations with round or elliptical bodies (Reference 1). During the initial 
development of program DEMON3 (Ref. 2), it was found necessary to model the finned 
section, comprising the lifting surfaces attached to the interference shell, with one set of 
constant u-velocity panels. As described in Reference 2, a convergent solution can not be 
obtained if the lifting surfaces in the finned section are modeled by constant u-velocity 
panels and if the interference shell in the same finned section is modeled by source or triplet 
panels. This fundamental problem appeared to be related to mathematical incompatibilities 
that prevent convergence in the matrix solution when the number of the two different types 
of panels is increased. This problem gave rise to Task 1 above. 

In the above task descriptions, program NOSEVTX refers to the vortex shedding 
computer program described in Reference 3. An improved version of program NOSEVTX, 
designated NOZVTX, is described in Reference 4. The final improved version of the vortex 
shedding program applicable to bodies with noncircular cross section is designated 
SUPVTX and is documented in Reference 5. The vortex shedding model employed for the 
body in DEMON3 makes use of the NOZVTX methodology in Reference 4 with some 
updates extracted from program SUPVTX of Reference 5. 
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During the development of the computer DEMON3 program, various portions of its 
emerging technology were extracted, modified, and/or improved resulting in various 
engineering level computer programs assembled to accomplish special tasks for NASA, U. 
S. Navy, and U. S. Army Agencies. A special version was developed for supersonic finned 
configurations with axisymmetric bodies. The NASA Langley Center version was 
designated LRCDM2 (Reference 6). The axisymmetric body version was extended to 
generate bulk force input for NASTRAN structural analysis resulting in the NWCDM- 
NSTRN program described in References 7 and 8. The NWCDM-NSTRN version can also 
include effects of steady angular rates in the calculated missile aerodynamic loadings. In 
addition, the LRCDM2 program was modified and the forebody vortex model improved to 
treat configurations with wrap-around fins in the AMICDM program (Reference 9). The 
DEMON3 program, without body vortical effects, was modified and extended to include 
supersonic inlets for the prediction of aerodynamic characteristics of airbreathing 
configurations at low angle of attack. This version was designated DM3INL and 
documented in Reference 10. The inlet modeling scheme, based on triplet panels, was 
extracted from DM3INL and improved which resulted in a separate module designated 
INLADD for estimating additive forces and moments associated with supersonic two- 
dimensional and axisymmetric inlets (Reference 11). Finally, as part of a DEMON3 
improvement effort, the SUBSAL program was developed for the Naval Weapons Center 
with NASA Langley Center as the monitoring agency. Program SUBSAL applies to 
subsonic configurations consisting of an axisymmetric body with up to two finned sections. 
This work is documented in Reference 12. 

The following sections describe the general approach embodied in program DEMON3, 
provide some details of the methods of analysis, and provide examples of the results 
predicted by the program. Limitations in the methodology of the DEMON3 program are 
discussed. This report includes a section describing the stepwise procedure incorporated in 
the DEMON3 computer program followed by descriptions of the input, output, and a 
sample case. Recommendations for improvements are furnished in the concluding section 
of this report. 
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LIST OF SYMBOLS 


b/2 exposed span of a lifting surface. See Figure 5 

c fin local chord 

c n span load, normal force (lbs) /unit span divided by dynamic head times local 

chord 

C i rolling moment coefficient about xg-axis, positive right fin down, looking 

forward, moment/(qJ> Ref L Ref ). See Figure 19 

C m pitching moment coefficient about yg-axis, nose up positive, 

moment/ (qjSRefLggf). See Figure 19 

C n yawing moment coefficient about Zg-axis, nose to right positive, 

moment/ (qjS Ref L Ref ). See Figure 19 

C N normal force coefficient along Zg-axis, force/ (q„S Re f). See Figure 19 

C Y side force coefficient along yg-axis, force/(qjS Ref ). See Figure 19 

C p pressure coefficient, (p-PoJ/q* 

D diameter of circular body or equivalent diameter of noncircular body based on 

cross sectional area 

FVN influence function, velocity component normal to surface induced by panel with 
unit strength 

L length of body 

L Re f reference length (normally body diameter) 

free stream Mach number 
p local static pressure 

p„ free stream static pressure 

q*, free stream dynamic pressure, 1/2 (p V„ 2 ) 

S Re f reference area (normally body cross sectional area) 

u,v,w axial, lateral, upward perturbation velocity components along body fixed 
• coordinates Xg, yg, Zg, refer to Figures 2 and 19 

V R magnitude of resultant velocity 

V OT free stream velocity 

x axial distance from body nose 



x B , yB' z B 


Xq P 

x mom 

X W'YW' Z W 


Yf 

a 

«c 

<t> 

r 

p 

<7 

Subscripts 


body fixed coordinate system with origin at nose tip, refer to Figures 2 and 19 

axial location of center of pressure 

moment center location measured aft of body nose 

wing coordinate system for a fin set parallel to Xg,y B ,z B system with origin on 
body longitudinal axis, xg, at leading edge of fin set 

distance from root chord along span of lifting surface. See Figure 5 

angle of pitch, degrees; impressed normal velocity 

included angle of attack, degrees, angle between free stream velocity vector and 
body longitudinal axis, xg, refer to Figure 19 

angle of roll, degrees, positive right wing down, refer to Figure 19 

vorticity; discrete vortex strength 

free stream density 

panel strength 


B 

BB 

F 

FF 

I 

II 
IF 
FI 
TE 


condition at infinity 
body 

body-on-body 

fin 

fin-on-fin 
interference shell 

interference shell-on-interference shell 
interference shell-on-fin 
fin-on-interference shell 
trailing edge 
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GENERAL APPROACH 


The need exists for fast, engineering level aerodynamic prediction methods capable of 
analyzing noneonventional wing-body-tail configurations. Methods based on experimental 
data bases are usually limited to geometrically restricted configurations; for example, 
conventional configurations with axisymmetric bodies and cruciform fin sections. Methods 
based on high order flow solvers (Computational Fluid Dynamics) contain the proper 
physics and can handle complex geometries but require substantial computer resources. 
Noneonventional configurations can be analyzed by panel methods since these methods are 
applicable to nonaxisymmetric bodies and arbitrary lifting surface layouts. Supersonic (and 
subsonic) panel methods are derived from the linear potential flow equation. The quality of 
the solution depends on the level of mathematical sophistication (treatment of boundary 
condition, treatment of continuity across the panels, etc.) encompassed in the particular 
panel method. High order panel methods provide the best solutions. However, high order 
panel methods can require as much computer resources as the simpler of the CFD-based 
methods, and the results will still be limited to the range of validity of linear theory. In 
supersonic flow, the best panel solutions can not inherently account for nonlinearities 
associated with the presence of shocks and vortical effects. 

The DEMON3 program described in this report makes use of low order supersonic 
panel solutions for the sake of economy. In order to increase the usefulness of the code, 
relatively simple models are incorporated to account for the nonlinear effects of lifting 
surface wakes and flow separation vortices from the the body which may be noncircular in 
cross section. The vortex tracking and vortex effects calculations are based on two- 
dimensional (slender body) concepts and are uncoupled from the three-dimensional 
paneling analysis. However, the vortical effects are included in the pressure calculations at 
points on the body and lifting surfaces. 

In addition, nonlinear effects of compressibility due to close proximity of shocks to the 
surfaces of the configuration can be optionally included. This is accomplished by calculating 
the pressures from the nonlinear two-dimensional shock expansion and Newtonian or 
impact pressure calculation methods. One of the unique features of DEMON3 is the 
extension of the two-dimensional nonlinear pressure calculations to include effects of 
aerodynamic interference obtained from linear theory (the panel solutions). 


METHODS OF ANALYSIS 


In order to compute pressure distributions, component loads, and overall forces and 
moments on a supersonic configuration consisting of a body which may have noncircular 
cross section and with up to two fin sets in various layouts, program DEMON3 makes use 
of different supersonic panel solutions to model the body and the fin sets. Specifically, the 
body is modeled by triplet panels developed by Woodward (Ref. 13) as an improvement to 
the supersonic source panels (Ref. 14). The lifting surfaces and the portion of body next to 
the lifting surfaces in a fin section are modeled by Woodward's constant u- velocity panels 
(Ref. 15) for lift on the surfaces and for lift carry-over on the interference shell, and thickness 
effects are modeled by supersonic planar source panels for thickness also originated by 
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Woodward (Ref. 14). A modified version of the vortex cloud method including the vortex 
tracking scheme described in References 4 and 5 is incorporated in DEMON3 to model 
vortex shedding from the body and to track vorticity down the length of the body. The 
concentrated multiple vortex wake model for the lifting surfaces as developed for the work 
reported in Reference 6 (the LRCDM2 computer program) is also incorporated into 
DEMON3. 

Program DEMON3 analyzes a given configuration from the nose to the base in a four 
(4) step procedure organized by an executive routine. The stepwise procedure is described 
in a later section describing the DEMONS computer program. In this process, body-on- 
lifting surfaces and lifting surfaces-on-body interference as well as vortical interference 
effects are included. For example, after the body solution is accomplished the effects of all 
the triplet panels are included in the flow tangency condition applied at points on the lifting 
surfaces. For each fin section, the strengths of the constant u- velocity panels on the lifting 
surfaces are obtained from a unified matrix solution which includes the constant u-velocity 
panels in the interference shell on the body next to the lifting surfaces. In this way, lift 
carry-over or lifting surface-on-body interference effects are modeled. Further details of the 
modeling schemes are given below. 


Lifting And Volume Panels For The Body 

In Reference 1 describing the predecessor (DEMON2) of DEMON3, the body is 
modeled by a distribution of supersonic source panels on its surface. The source panel 
models both lift and body volume effects. The particular source panel solution was 
developed by Woodward and is described in Reference 14. This type of low order body 
surface paneling generates a multitude of two-dimensional disturbances radiating towards 
the interior of the body. In the application to long bodies with open bases or to complete 
closed bodies, the wave-like disturbances reflect back and forth and can cause numerical 
difficulties in the solution of the source panel strengths especially at low supersonic Mach 
numbers. As a result, the magnitudes of the resulting panel strengths can be erratic and 
very high thereby causing the perturbation velocities to be unrealistic at points on and off 
the body although the flow tangency boundary condition is still satisfied. The distribution 
of surface pressures based on the velocity components calculated at points on the body can 
also be erratic and unrealistic. In an effort to remedy this difficulty. Woodward developed 
the triplet panel which essentially eliminates the two-dimensional portion of the internal 
waves. The salient features of the triplet panel are summarized next. Details of the solution 
are given in References 13 and 16. 

Triplet Panel Singularity 

The triplet panel singularity is a linear combination of source and vortex singularity 
distributions. Since the vortex (or doublet) sheet may be considered the result of combining, 
in a special limiting process, two source sheets of equal and opposite strength, the 
combination of a vortex sheet and a source sheet has been termed a triplet. This 
superposition is illustrated in Figure 1(a) which shows in the two-dimensional sense that the 
axial and normal velocity components add in the flow field above the panel, but cancel 
exactly below the panel, resulting in the desired unidirectional perturbation velocity field. 
The triplet panel models lift and body volume effects. 
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As described in References 13 and 16, the extension of the above two-dimensional 
unidirectional concept to a body having arbitrary cross section requires a triplet panel 
grouping of six (6) panels in which the vortex loops are closed. A typical example is 
indicated in Figure 1(b). Depending on the orientation, the panels within a grouping are 
designated axial, circumferential, and radial triplets. The references show various pressure 
distributions calculated with the new triplet panel approach and by the original source 
panel solutions. In most cases, the erratic behavior of the original source distribution 
especially near the base of the bodies has been eliminated. 

In the application to bodies, the triplet panel solution originally set up by Woodward 
made use of the simplifications associated with the assumption of a vertical plane of 
symmetry with regard to both geometry and the oncoming flow. In the DEMON3 program, 
the symmetry built into the triplet solution with regard to the oncoming flow had to be 
relaxed in order to handle nonzero roll angle. This was accomplished by modifying the 
subroutines of a later version of Woodward's USSAERO code which includes the triplet 
panel update. The original version of the USSAERO code for subsonic and supersonic 
analysis of aircraft configurations is described in Reference 17. 

Flow Tangency Condition On The Body 

The strengths of the triplet panels used to model the body of a complete configuration 
are obtained from the set of simultaneous equations that are the result of applying the flow 
tangency condition at the control points of the triplet panel. The control points are located 
at the centroids of the individual triplet panel areas. One control point is shown in Figure 
1(b). At a given control point, the velocity component normal to the plane of one triplet 
panel induced or influenced by another triplet panel of unit strength is expressed in terms of 
the influence function FVN. As mentioned earlier, the solution and thus the influence per 
unit strength of a triplet panel involves groupings of six (6) panels. In matrix notation, the 
flow tangency condition can be written as 

[FVN BB ] [<t b ] = [ag] (1) 

where subscript BB denotes body-on-body triplet panel influence, each row in the matrix 
FVN represents the sum of the influences per unit strength of all the triplet panels at the 
control point of one panel, vector [tr B ] contains the strengths of the panels, and vector [a B ] 
contains the set of impressed velocities at the control points. 

In program DEMON3, the impressed normal velocities depend on the free stream and 
the panel inclination angle. 


Lifting And Thickness Panels For The Fin Section 

A typical fin section made up of cruciform fins on an axisymmetric body is shown in 
Figure 2. As stated earlier, program DEMON3 is not limited to this conventional layout and 
the following discussion pertains to nonconventional layouts with nonaxisymmetric body 
cross sections as well. 
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The lifting surfaces and the interference shell indicated in Figure 2 are covered by 
constant u-velocity panels originally called constant pressure panels by Woodward and 
described in Reference 14. These panels model the lift effects and are positioned in the 
chordal planes of the lifting surfaces. Therefore, the constant u-velocity panel method is not 
a surface paneling method. 

In addition, thickness effects of the lifting surfaces are modeled by a layout of 
supersonic planar source panels. In the DEMON3 program, the spanwise layout of the 
source panels is the same as for the lifting constant u-velocity panels. However, in order to 
resolve streamwise thickness effects, the chordwise layout is normally more dense. The 
strength of a source panel is directly related to the local thickness slope of the lifting surface. 

Constant U- Velocity And Source Singularities 

The particular singularity solution of the constant u-velocity panel is one for generating 
lift in supersonic flow. This means that the axial and lateral velocities immediately above 
and below the panels are equal and opposite in sign. The normal component induced by the 
panel is continuous across the plane of the panel and is made to counteract the impressed 
normal velocity. 

The source panel singularity has continuous axial and lateral velocity components 
immediately above and below its plane. The normal velocity component is continuous 
across the plane of the source panel. The source panel strengths are directly related to the 
local thickness slope. 

A detailed description of the two types of panel solutions and the particular solutions 
in the respective panel planes are listed in Appendix II of Reference 18. 

Flow Tangency Condition On The Lifting Surfaces 

The flow tangency condition is applied at the control points of the constant u-velocity 
panels on the lifting surfaces and on the interference shell. Two control points are indicated 
in Figure 2. 

On the lifting surfaces, the impressed velocity in the flow tangency condition includes a 
contribution from the free stream, body-induced flow (upwash), and effects of external 
vortices from the upstream fin section (if applicable) and the body. The constant u-velocity 
panels on the lifting surfaces interact with the constant u- velocity panels on the interference 
shell. In this way, lift carry-over effects are included. On the interference shell around the 
body, the impressed velocity is due to lifting surface thickness only. All other impressed 
velocities are included in the triplet panel solution for the body as described in an earlier 
section. 

The velocity component induced by one constant u-velocity panel of unit strength 
normal to the plane of another panel at its control point is the influence function FVN. The 
boundary condition can be expressed in matrix notation as follows: 
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( 2 ) 


FVNpp FVN 1f ' 


V F ‘ 


a F 

FVNpi FVN n 




_ a I . 


In the above expression, subscript FF denotes one constant u-velocity panel on a lifting 
surface influencing another panel on the same or other lifting surface, subscript II refers to 
one panel on the interference shell influencing another panel on the interference shell, 
subscript FI means one panel on a lifting surface influencing a panel on the interference 
shell, and subscript IF relates to one panel on the interference shell influencing a panel on a 
lifting surface. Vector [a] contains the strengths of the constant u-velocity panels on the 
lifting surfaces and the interference shell in a fin section. Vector [a] represents the set of 
impressed velocities enumerated above. Note that [aj] is zero unless the thickness effects of 
the lifting surfaces are included. 


Body Vorticity Model 

The body vorticity model incorporated in program DEMON3 is nearly identical to the 
engineering level vortex shedding prediction method described in Reference 4. One of the 
differences is the replacement of the source panels by triplet panels for the inviscid 
modeling of the body. In addition, information including vortex strengths and positions as 
well as body forces and moments are exchanged between the body vorticity routines and 
the fin section routines. 

In short, the body is represented by the supersonic triplet panel method (as described 
in an earlier section), and the leeside vortex wake is modeled by discrete vortices in cross 
flow planes. The three-dimensional steady flow problem is reduced to a two-dimensional, 
unsteady, separated flow problem for solution which makes use of conformal mapping 
schemes to transform noncircular body cross sections to a circular contour. The predicted 
pressure distribution on the body under the influence of free stream and the body 
separation wake is used to calculate the body loads. The effects of the body separation 
wake are included in the analysis of the fin section. A typical example of the vortex patterns 
on a 3:1 elliptic body is shown in Figure 3, taken from Reference 4, for zero and 45 deg roll 
angles, respectively. 


Bernoulli Pressures, Forces And Moments 

The calculations of pressure at points on the body and at points on the fin sets are 
based on the compressible Bernoulli pressure/ velocity relationship. Due to the presence of 
vortices on the body when the angle of attack is sufficiently high, the formulation of the 
pressure equation for the body is somewhat different than that employed on the fin set. 
Details of either pressure calculation are available in Reference 4 and in Reference 6 for the 
body and for the fin sets, respectively. For the sake of completeness, the pressure 
calculation methods are summarized below. 
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Pressures On The Body 

The following pertains to the section(s) of the body which are not part of the fin sets of 
a complete wing-body-tail configuration. Circumferential pressure distributions are 
required at many axial stations along the body in order to calculate the contributions to the 
overall aerodynamic forces and moment acting on the complete configuration. The pressure 
distributions on the body sections also determine the flow separation points in the cross 
flow planes. At each axial station, the pressures are calculated at 73 points placed at equal 
angular intervals on the body contour. The Bernoulli equation is written in the following 
pressure coefficient form where y = 1.4 for air. 


C 


P 


P - P„ 


% pv 




y- 1 
2 



4 


(3) 


The incompressible contribution C Pl 


is given by the following expression. 



1 



2cosa 

c dip 
V dx 


(4) 


In the above, U is the total velocity including free stream at a point on the body. In the 
calculation of this total velocity, the axial and lateral perturbation velocity contributions 
from the triplet panels (modeling body lift and volume effects) are included. The lateral 
velocity contributions from the external vortices, composed of the wake vortices from the 
upstream fin set and the body shed vortices, are added to the body-induced velocities. The 
last term in Equation 4 represents the axial velocity component induced by the external 
vortices by virtue of their inclination relative to the body centerline. The details of this 
contribution to the axial flow component are given in Reference 4 which shows that this 
component is related to the coordinates of the vortices in successive cross flow planes. Since 
the pressure coefficient as expressed in the above formulations depends on velocity 
components obtained from linear theory, it is possible that local pressures lower than free 
stream static pressure are calculated. Therefore, the calculated pressure coefficients are 
limited to the free stream value as follows (y = 1.4 for air). 




(5) 


It is worth noting here that the vortex tracking and vortex-induced effects calculations 
are based on slender body theory (for example, two-dimensional doublets). This means that 
apart from the separation point determination procedure, the three-dimensional body flow 
model (the triplet panel solution for the body at angle of attack) and the fin set flow model 
(the constant u- velocity panel and thickness panel solutions) do not take part in the vortex 
tracking calculation. The exception is the axial velocity component from the body volume 
effects (zero angle of attack triplet solution) which is included in the vortex tracking 
procedure. In that procedure, the vortices are made to follow the streamlines associated 
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with the flow which includes all upstream effects. In this way, the vortex model is 
essentially uncoupled from the panel method modeling scheme for the body. On the fin set, 
the effects of the vortices are included in the flow tangency boundary condition as described 
below. 

The individual aerodynamic force coefficients are obtained by multiplying the 
pressures calculated by the method above times the appropriate areas on the body surface. 
The individual moment coefficients are calculated by multiplying the individual force 
coefficients by the appropriate moment arms. The total body contributions to the overall 
force and moment coefficients are then obtained from the sums of the individual force and 
moment coefficients. 


Pressures On the Fin Set 


The pressures are calculated immediately above and below the plane of the constant u- 
velocity panels at the individual panel area centroids. These panels are distributed on the 
lifting surfaces and on the interference shell in the fin set. On the interference shell, the 
pressures are calculated on the exterior of the panels only. 

The Bernoulli pressure /velocity relationship used in the fin set is given by the 
following expression in terms of the resultant velocity V R (y = 1.4 for air). 






oo 


1 



1 


The resultant velocity ratio is given by 


( 6 ) 


V 2 2 2 

R „.2u 2v, , A , 2w . a , u + v + w 

— x “ 1 + — cosa - — smct sin# + — sina cos# + = (7) 

„2 V c V c T V c „2 

V oo oo oo V 

oo oo 

where a c is the included angle of attack and # is the angle of roll related to the the pitch-roll 
definition for angle of pitch and side slip (refer to Ref. 6, Equation (3) ). The perturbation 
velocities u, v, w induced by the body, lifting surface, and interference panel methods can 
be unduly large in magnitude and cause the term in the square brackets in Equation (6) to 
become negative. In such cases, the pressure coefficient is limited to the value given by 
Equation 5 shown above. In program DEMON3, therefore, the pressures calculated from 
the panel solutions with the Bernoulli expressions are nonlinear in terms of angle of attack 
by virtue of its formulation and the imposed limiting value. 

The perturbation velocities u, v, and w included in the pressure calculation at points in 
the fin set (Eqs. 6 and 7) are made up of contributions from the body triplet panels, the 
constant u-velocity panels and thickness panels on the lifting surfaces, the interference shell 
of the fin set, and from the vortices generated by the upstream body section and fin set if 
applicable. In this procedure, the vortices are "frozen" at the leading edge of the fin set 
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(actually the axial station corresponding to the leading edge of the fin rootchords). This 
means that at the leading edge position, the vortices are considered as two-dimensional 
point vortices with known strengths and positions calculated by the tracking scheme. This 
is equivalent to assuming the vortex paths to be parallel to the body centerline through the 
fin set section. This is in contrast to the vortex tracking scheme in the DEMON2 program of 
Reference 1 in which the vortices are tracked through the fin set section. It was found in 
later developments that the DEMON2 tracking approach adds complications without 
necessarily increasing accuracy. 

Special care has to be taken in the calculation of the lifting surface panel contributions 
to the perturbation velocities used in the pressure calculations immediately above and 
below the surfaces of the lifting surfaces. Because of the planar (non-surface panel) nature 
of both the constant u-velocity and the thickness (or source panel solutions), the DEMON3 
program has to determine the proper values of the panel-induced velocities immediately 
above and below the panel planes. 

In the fin set pressure calculations, the external vortices (assumed to run parallel to the 
body centerline through the fin set section) only induce lateral (v, w) velocity components. 
It should be noted, however, that the lateral velocities induced by the vortices are included 
in the flow tangency condition applied at the control points of the constant u-velocity panels 
distributed on the lifting surfaces as described in an earlier section. Therefore, the constant 
u-velocity panels under the influence of the vortices will induce u, v, w velocity components 
which reflect the presence of vortices to first approximation. This approximate approach is 
successfully employed in the various derivations of the DEMON3 program (Refs. 6-9, 12) as 
well as in other aerodynamic interference prediction methods based on the equivalent angle 
of attack approach. 

The force and moment calculations are performed in essentially the same manner as 
used on the body with the exception that the areas over which the pressures act are the 
actual constant u-panel areas. 


Shock Expansion And Newtonian Pressure Calculations 

As an option, program DEMON3 can calculate pressures on both the body and the 
lifting surfaces based on two-dimensional shock expansion and Newtonian theories. The 
designation two-dimensional means that these nonlinear pressure calculations are 
performed along meridians of the body and along chord wise strips on the lifting surfaces. 
The optional pressure calculation methods are applicable to configurations at high Mach 
numbers (in excess of approximately Mach 2.5 up to Mach 6). The Newtonian method is 
implemented in DEMON3 in its simplest form and assumes zero pressure for any "shaded" 
portion of the configuration. It should only be applied to configurations for which the Mach 
number is in excess of 5. 

The shock expansion calculations are based on the assumption of attached shocks, and 
the strip theory approach used in DEMON3 is general as opposed to second order. In the 
general shock expansion theory, the pressure is constant along a segment of a strip. In 
second order theory, the pressure varies along the segments of the strip. The validity of 
general shock expansion versus second order depends on the hypersonic similarity 
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parameter. This parameter is the ratio of free stream Mach number over fineness ratio. If 
this ratio is much larger than one (high Mach number), general shock expansion theory is 
valid. If the ratio is about one, second order theory should be used. 

In programs such as DEMON3 and its derivatives based on either low order panel or 
line singularities to model the body, linear theory requires that the Mach cone attached to 
the body nose lie outside the body meridional contour as modeled by the respective 
programs. This requirement can be violated in any case when the body nose is blunt 
and/ or when the free stream Mach number is high for all but the sharpest body noses. This 
problem can be circumvented by reducing the body slopes near the nose. When this user 
arranged procedure is used, the pressures (Bernoulli, as well as the shock expansion and 
Newtonian) near the nose will be underpredicted, however. 

The unique feature of the shock expansion and Newtonian pressure calculation 
methods in DEMON3 is the inclusion of an approximate scheme to correct the nonlinear 
pressures for strip-on-strip, fin-on-fin, mutual body-fin, and vortical aerodynamic 
interference effects. In fact, the interference on the flow deflection angles is obtained from 
the panel solutions (linear theory) and the vortical induced effects. The modified flow 
deflection angles are used in the nonlinear shock expansion and Newtonian pressure 
coefficient calculations. This engineering level procedure was originated by Carlson (Ref. 
18) and first implemented in the LRCDM2 program (Ref. 6). The actual formulations of the 
nonlinear pressure calculation methods and the corrections to the flow deflection angles are 
described in detail in the main body and in Appendix D of Reference 6. A summarized 
account is provided in Reference 19. In the LRCDM2 program, the methodology is really 
applicable to the lifting surfaces only because analytical tangent wedge (oblique shock) 
solutions are incorporated for the pressures and local Mach number immediately behind the 
shock on the leading edge. The LRCDM2 program can also calculate shock expansion and 
Newtonian pressure as an option on the axisymmetric forebody. However, the analytical 
tangent wedge (oblique shock) expressions are used to compute pressure at the nose tip as a 
first approximation instead of using the proper tangent cone values which are not available 
in exact analytical form. This limitation has been relaxed somewhat in the DEMON3 
program by the inclusion of approximate tangent cone relationships. The procedure and 
specific shock cone relationships used in program DEMON3 are enumerated below. 

Tangent Cone Relationships Used In DEMON3 

In program DEMON3, the optional shock expansion and Newtonian pressure 
calculation method including the flow deflection angle corrections are applicable to the 
forebody which may be noncircular in cross section. The geometric layout scheme in 
DEMON3 establishes the meridional strips and the axial stations from the body nose to the 
forward fin set at which the optional pressure calculations will be performed. At the 
beginning or leading edge of each strip (at the body nose), the analysis assumes the strip 
under consideration to be part of an axisymmetric cone and independent from the other 
strips around the body. For each strip at the body nose, the shock expansion analysis 
requires values of shock detachment angle (in order to determine if the analysis can 
proceed), Mach number normal to the shock, cone surface Mach number, cone surface 
pressure, and ratio of stagnation pressure across the shock. These parameters are not 
available in exact analytical form as is the case for the tangent wedge case. One possibility is 
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to incorporate data bases or to employ approximate relationships. The present version of 
DEMON3 uses various expressions that approximate the tangent cone shock parameters. 

The shock detachment angle and the cone surface Mach number immediately behind 
the shock are determined using the approximate solutions developed by Hammitt and 
Murthy and described in Reference 20. These approximate solutions are based on Taylor 
series expansions for the velocities between the cone shock and body and involve the body 
cone half angle, free stream Mach number, and shock angle. The value of the Mach number 
normal to the shock and the cone surface pressure are obtained from the tangent cone 
approximations modified by Edwards and specified in the appendix of Reference 21. These 
approximations are based on previously developed empirical representations and involve 
the flow deflection angle and free stream Mach number. In program DEMON3, the 
approximate tangent cone relationships are programmed in subroutine SHKXB. 


Lifting Surface Wake Vortex Model 

In the DEMON3 program, the wake at the trailing edges of the lifting surfaces or fins in 
a fin set is represented by a set of discrete vortices. In addition to the trailing edge vorticity 
related to attached flow, the wake model includes leading and side edge vortices which 
develop and dominate the fin wakes as the angle of attack is increased. For lifting surfaces 
with long side edges, vorticity can be generated along the edge for angles of attack as low as 
5 deg. A typical fin wake model as generated by the program is depicted in Figure 4 which 
shows the trailing edge vortex (more than one can exist) and the combined leading and side 
edge vortex which is elevated from the fin plane, These vortices are included in the 
tracking scheme that is part of the body vortex shedding model described in an earlier 
section. Therefore, the fin wake vortices influence the body vortex shedding process along 
the afterbody as well as the loads acting on the tail fin set. 

The actual wake model used in DEMON3 is nearly identical to that of the LRCDM2 
program (Ref. 6), and a summarized description is also available in Reference 19. The only 
difference in the wake model is related to the attached flow vortices as described below. 

Updated Trailing Vortex Model For Attached Flow 

As described in References 6 and 19, the trailing edge vorticity r T g associated with 
attached flow on the lifting surface can be related to the span loading c n as follows (c is local 
chord). 


1 _ 

V 



1 _a_ 

2 ay. 


<cc n ) 


( 8 ) 


The dted references contain the analytical expressions derived from the above differential 
equation on the basis of only one extremum in the span load distribution. An example is 
shown in Figure 5 for a delta fin on a body and under the influence of an external vortex. In 
actuality, the span loading may exhibit multiple extrema depending on the oncoming flow 
and possibly the particular fin-on-body attachment when the body is noncircular. 
Therefore, the integration of Equation (8) was extended to account for multiple extrema. 
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The strengths and positions of the wake vortices are calculated in subroutine SPNLD of 
program DEMON3. 


RESULTS 


This section presents results calculated with the DEMON3 aerodynamic prediction 
program. Results have been obtained for a variety of configurations. Comparisons are 
made between DEMON3 results, experimental results, and results from other prediction 
programs and theory for the purpose of validation and verification. The first comparison is 
for a conventional missile configuration with a round body and cruciform canard and tail 
fin sections (Ref. 22). The second and third configurations are monoplane wing and tail 
configurations with noncircular bodies (Refs. 23 and 24). Finally, the shock expansion 
option in the code is exercised and compared with data for a rectangular wing (Ref. 25) and 
for an ogive cylinder body (Ref. 26). 

Cruciform Canard Control, Cruciform Tail Configuration 

Comparisons of measured and predicted characteristics of a canard-body-tail model 
are described in this section. The configuration is the canard-controlled wind tunnel model 
described in Reference 22 and depicted at the top of Figure 6. The body is a 2.25-caliber 
ogive nose followed by a 21.9-diameter cylindrical body. The canard fins have an aspect 
ratio of 2.00, a taper ratio of 0.30, and a body radius to fin semispan ratio of 0.226. The 
canard leading edge is located 3.73 diameters aft of the nose tip. The tail fins under 
consideration have an aspect ratio of 1.06, a taper ratio of 0.58, and a body radius to fin 
semispan ratio of 0.19. The tail fins are located 12.7 diameters aft of the canard trailing edge. 
Overall loads on the model at several supersonic Mach numbers are available, but only 
results for M^ = 2.5 are presented here. 

Measured and predicted longitudinal aerodynamic characteristics are indicated in 
Figure 6 for the canard-body-tail model shown at the top of the figure at zero roll angle. 
Figure 7 shows the measured and predicted lateral directional aerodynamic characteristics. 
The Mach number is 2.5. The right and left horizontal canard fins are deflected +5 deg and 
-5 deg, respectively, for right wing up (negative) roll control. In addition to the DEMON3 
predictions, results are included from the SUPDL program (designated NWCDM-NSTRN in 
Ref. 8). Program DEMON3 models the body of the configuration with supersonic triplet 
panels, whereas program SUPDL makes use of linearly varying supersonic line 
sources/sinks and doublets. Both programs employ the same fin set modeling method, and 
both include models to account for the nonlinear effects of body vorticity and fin wakes as 
described in earlier sections. In Figures 6 and 7, SUPDL results are shown with 
(w/BDYSHD) and without vortex shedding from the length of body between the canard 
and tail fins. All DEMON3 results include body vortex shedding. 

The normal force and pitching moment coefficients shown in Figure 6 exhibit some 
nonlinearity with angle of attack. All three predictions shown in Figure 6 show the same 
nonlinear trend due to the inclusion of the fin wake effects. Body vortex shedding has 
negligible influence in the longitudinal charateristics. In general, the DEMON3 results agree 
best with the experimental data. 
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The yawing moment, side force, and rolling moment coefficients shown in Figure 7 are 
small compared with the longitudinal characteristics. The rolling moment behavior is of 
major interest, however. Experimental data for the rolling moment acting on the 
configuration with the tail section removed are also shown in the bottom portion of the 
figure. For the complete configuration at low angles of attack, the canard rolling moment 
(tail off) is counteracted by the rolling moment induced by the canard wake on the tail fins 
resulting in near zero overall rolling moment. As the angle of attack is increased, this 
adverse effect disappears. Above about 10 deg angle of attack, the magnitude of the overall 
rolling moment is higher than the canard fin alone (or tail off) value. At angles of attack in 
excess of 15 deg, the overall rolling moment returns to the canard alone value up to canard 
fin stall condition resulting in zero rolling moment. The DEMON3 and SUPDL predictions 
match the tail off value. The nonlinear trends of the overall rolling moment are predicted 
fairly well by the DEMON3 and SUPDL codes provided body shedding is included in the 
latter. The lines associated with the predictions connect points calculated at 5 deg intervals 
in angle of attack. 


Mid Monoplane Wing and Tail Configuration 

This section presents comparisons for the monoplane wing and tail model of 
Reference 23 depicted in Figure 8. The body of the configuration is nearly circular over the 
ogive nose and transitions to an elongated (vertical) circular section on the forebody. Before 
the wing leading edge, the cross section transitions to a circular bottom and square top. 
Over the wing root chord, the cross section transitions back to an elongated (vertical) circle 
and maintains this shape to the base. Tails are located at the rear of the configuration. 

Comparison between measured and predicted results for = 1.6 is shown in Figure 9. 
Results are also included from preliminary analysis and design programs MISL3 (Ref. 27) 
and Missile DATCOM (Ref. 28). For MISL3 and Missile DATCOM, a circular body has been 
used. The normal force coefficient is predicted well by all the programs. Figure 9(a). 
Program DEMON3 predicts the pitching moment well for angles of attack up to ten degrees. 
Figure 9(b). Above ten degrees the pitching moment may be in error because of wing 
and/or tail fin stall. DEMON3 does not currently include a stall model. The MISL3 and 
Missile DATCOM programs do not predict the pitching moment indicated by the 
experimental data. Note that for this configuration, the pitching moment from the nose and 
tail fins counteract one another; therefore, the pitching moment is sensitive to the wing 
center of pressure. The wing is large but located near the moment center; therefore, an error 
in the wing center of pressure can have a large effect on the pitching moment. 

There are several possible reasons why MISL3 and Missile DATCOM do not predict 
the correct pitching moment: the wing center of pressure is not predicted correctly and/or 
the effect of the noncircular body affects the pitching moment (not modeled by MISL3 and 
Missile DATCOM). The body cross section grows from an elongated (vertical) circle to a 
circular bottom with a square top cross section before the wing leading edge. Over the wing 
root chord section it transitions back. This transition back to the elongated (vertical) circle 
accelerates the flow over the top portion of the body and wing. This accelerated flow affects 
the wing load, load distibution, and hence pitching moment. MISL3 and Missile DATCOM 
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do not include this effect. The triplet panel body model in DEMON3 models the effects of 
the actual body shape. 

The measured and predicted locations of the configuration center of pressure 
(measured from the moment center and referenced to the mean aerodynamic chord of the 
wing) are shown in Figure 9(c). DEMON3 predicts the center of pressure location well for 
this configuration. 

Comparison between measured and predicted results for M„ = 2.0 is shown in 
Figure 10. The DEMON3 results overpredict the normal force and pitching moment at the 
higher angles of attack for this Mach number. This may be due to wing and/ or tail fin stall 
which is currently not modeled by the program. Less wing and/ or tail fin normal force 
would reduce both the overall normal force and pitching moment. 

High Monoplane Wing Tri-tail Configuration 

Comparisons of DEMON3 with predicted results from a modified version of the MISL3 
program for the Pegasus™ launch vehicle (Ref. 24) are described in this section. The 
configuration consists of a circular cross section forebody which transitions to noncircular 
cross section ahead of the high wing. The noncircular body shape extends aft of the wing 
trailing edge. The body transitions back to circular before the tri-tail configuration. There 
are two horizontal tail fins and one vertical fin. This configuration is depicted in Figure 11. 

Overall pitch plane loads on the model at = 2.0 are shown in Figure 12 with and 
without 10 deg of horizontal tail deflection. The DEMON3 predicted normal force and 
pitching moment coefficients agree well with those obtained with a modified version of the 
MISL3 program. The effect of deflection also agrees well. One of the advantages of the 
DEMON3 program is its capability of providing aerodynamic force distributions which can 
be incorporated in a structural analysis. 


Shock Expansion and Newtonian Theory Applications 

In this section, two examples are discussed of shock expansion and Newtonian 
pressure calculations including interference corrections. The first example involves a 
rectangular wing alone, the second example is concerned with an ogive cylinder body. 

Rectangular Wing 

Pressure distributions are available (Ref. 25) for the beveled, aspect ratio 2, rectangular 
wing shown in Figure 13. The figure shows the pressure tap locations and the panel layout 
on one-half of the wing for the calculations referred to below. The chordwise pressures 
measured on the upper and lower surfaces at the one-half semispan location of the wing are 
shown in Figure 14. The angle of attack is 10.3 deg and the Mach number is 2.86. The figure 
shows the strong effects of the beveled portions of the wing as indicated by the 
experimental data. 

Comparisons between results calculated by the LRCDM2 program and the 
experimental data are shown in References 6 and 19. The purpose of the comparisons 
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shown in Figure 14 between the DEMON3 results and the experimental data is to ascertain 
that the optional nonlinear pressure calculation methods perform properly in the DEMON3 
program. The methods for the wing pressures in DEMON3 are essentially the same as in 
LRCDM2 with the exception of an improvement in the determination of the panel-induced 
velocities immediately above and below the plane of the panels. It is worth noting that the 
LRCDM2 program can handle a wing alone. The present version of the DEMON3 program 
requires the presence of a small body in "wing alone" calculations. 

The linear theory results generated by the DEMON3 program are based on a 10 
chordwise by 10 spanwise layout of constant u-velocity panels to model lift, and the same 
layout of source panels used to model the thickness envelope. The shock expansion and 
Newtonian analyses make use of 10 chordwise strips with 10 segments in each. A length of 
body with radius equal to the wing thickness is placed on the root chord. The results 
labeled "linear" are based on the panel solutions (linear theory) and the compressible 
Bernoulli pressure expression (Eq. 6). The results labeled "shock expansion, corrected" are 
based on the shock expansion theory with flow deflection angles corrected for strip-on-strip 
interference by the panel solutions as summarized in an earlier section. The results labeled 
"Newtonian, corrected" are based on the simple impact theory method (Cp = 2 sin 2 8 , 
Cp = 0.0 for shaded areas) with the flow deflection angle 8 corrected by the panel solutions. 

Because of the presence of the strong oblique shock attached to the leading edge, the 
linear result underestimates the pressure coefficients on the upper and lower surfaces of the 
wing. The corrected shock expansion results are in good agreement in the region near the 
leading edge. This case is representative of nonlinear compressibility influencing the 
pressures due to the close proximity of the leading edge shock to the wing surfaces. Over 
the flat portion, the linear and shock expansion methods agree well with the experimental 
data. On the beveled trailing edge, the linear and corrected shock expansion method 
predict lower than measured pressure on the upper (suction) surface. As mentioned in 
References 6 and 19, this is most likely due to boundary layer separation effects. On the 
lower (windward) side, the corrected shock expansion method matches the experimental 
data better. The corrected Newtonian result severely underestimates the pressures on the 
lower surface, and the pressures are zero on the upper surface. The Mach number in this 
case is too low for application of the Newtonian pressure method. 

In this example, the correction to the flow deflection angles is minimal due to the two- 
dimensional flow character in the midspan region of the rectangular wing. The basic strip 
theory is adequate already. The DEMON3 results shown in Figure 14 closely match the 
results generated by LRCDM2 published in References 6 and 19. 

Ogive Cylinder 

In order to assess the capability of the pressure calculation methods implemented in 
the DEMON3 program in the application to bodies, comparisons with experimental data for 
an ogive cylinder body are shown in Figure 15. The geometric details of the body are 
shown in the lower portion of the figure. The angle of attack is zero and the Mach number 
is 2.96. 

On Figure 15, experimental pressure coefficients are shown by the open symbols for the 
upper and lower meridians of the pointed ogive-cylinder. The two measurements should 
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be identical for the zero angle of attack case. The pressures overshoot the zero level in the 
region of the tangency point. 

In the application to a body alone, program DEMON3 does not need to model a 
dummy wing unlike the LRCDM2 program (Ref. 6, Section 3.5). The DEMON3 predictions 
are based on a distribution of triplet panels on the body consisting of 15 rings along the 
body with 8 panels each around the circumference. The rings are clustered near the body 
nose for the sake of resolving details of the nose shape. For the shock expansion and 
Newtonian pressure calculation methods, program DEMON3 automatically lays out 73 
strips around the body, and the program input was arranged to specify 13 axial stations at 
which circumferential pressure distributions are calculated. 

In Figure 15, the Bernoulli pressures calculated by Equation (3) with velocities induced 
by the triplet panels underpredict the experimental pressures from the nose tip up to about 
half the length of the ogive nose. Aft of this station, the Bernoulli pressures agree fairly well 
with the experimental data including the overshoot. The uncorrected shock expansion 
pressures match the experimental data considerably better on the nose although the 
pressures are somewhat underpredicted in this region. The corrected shock expansion 
pressures are a little lower still. Aft of the tangency point, the shock expansion pressures 
systematically underpredict. This deficiency is thought to be related to the use of general 
shock expansion rather than second order shock expansion. In general, on the nose these 
results are a definite improvement over the results obtained by the LRCDM2 code as 
published in Reference 6. The main reason for the improvement is the incorporation of the 
approximate tangent cone relationships in lieu of the tangent wedge expressions for 
conditions immediately behind the shock. The Newtonian prediction underestimates the 
pressures over the entire nose. In this case, the Newtonian-based pressures are not 
applicable due to the low Mach number. Additional comparisons with experimental data 
for noncircular bodies at angle of attack need to be performed. 


LIMITATIONS IN THE METHODOLOGY 


The DEMON3 computer program is subject to certain limitations and deficiencies due 
in part to the low order supersonic panel modeling employed and the associated limitations 
in the geometric details that can be modeled. Additionally, there is a shortcoming due to 
the general shock expansion theory employed in the optional nonlinear pressure 
calculations on the body. Five limitations and deficiencies are listed below. 


Supersonic Paneling Deficiencies 

The following deficiencies are related to the triplet panel and constant u-velocity panels 
used to model the body and the lifting surfaces, respectively. 
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Body Triplet Panels 

As stated in the technical approach, the DEMON3 program makes use of Woodward's 
triplet panels (Ref. 13) to model the body which may be noncircular in cross section. There 
are at least three problems associated with the triplet panel in its present form. 

(1) The triplet panel geometrical layout and theoretical solution are based on 
the panel being flat so that the four panel comers lie in the same plane. In 
cases for which the body cross sectional contours are similar from the nose 
to the base, all of the triplet panels will be planar and the solution will be 
valid within the constraints of linear supersonic panel theory and the 
specific characteristics of the low order triplet panel. If the body contours 
are dissimilar, the triplet panels may be nonplanar or twisted if they are 
defined with four comers. The triplet routines in DEMON3 (taken from a 
later version of the USSAERO program described in Reference 17 and 
modified for DEMON3 to handle asymmetric flow conditions ) will still 
generate solutions. It is clear that the solutions in cases with twisted panels 
may not be valid. However, a twisted panel can be divided into two 
triangular planar panels. Unfortunately, the present triplet panel layout 
routines do not divide four cornered twisted panels into three cornered 
panels automatically. However, the input of DEMON3 allows the user to 
arrange for two corners to coincide resulting in a triangular panel. 
Therefore, to the extent possible the user should avoid generating layouts 
resulting in twisted or nonplanar panels This can be done by rearranging 
and/or dividing a nonplanar panel into two triangular panels by judicious 
selection of panel layout coordinates in the input. 

(2) Notwithstanding the improvements embodied in the triplet panel relative to 
the source panel solution to avoid internal reflections (Refs. 13 and 16), it is 
possible to encounter cases where the panel solution results in erratic and 
very high panel strengths towards the base of the body. The flow tangency 
boundary condition is still satisfied but the panel strengths have large 
magnitudes which may oscillate down the body. As a result, the velocity 
components assume large values which affect the pressures calculated on 
the body. This problem can happen when the body is slender (or long) and 
the Mach number relatively low. In such cases, the residual or uncancelled 
portion of the influences radiating aft inside the body from the comers of 
the triplet panels along Mach waves or cones are amplified as they reflect 
back and forth. Other than updating to a higher order paneling method 
which mathematically assures that there are no internal influences, there is 
no easy remedy for the residual internal reflection problem related to the 
triplet panel. 

(3) The triplet singularity solution is based on the requirement that the triplet 
panel must lie inside the Mach wedge and Mach cones attached to the 
leading edge and its corner points. The wedge and cone angles depend on 
the free stream Mach number. At and/or near the nose, the panel 
inclination angles may be too large depending on the body nose geometry 
resulting in superinclined panels. In such cases, the DEMON3 program will 
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stop with a message described in a later section concerned with the 
computer program description. The user should reduce the panel 
inclination angles by lengthening the first sets of circumferential panels 
until the planes of the panels lie inside the Mach cone attached to the nose. 

Constant U-velocity Panels on the Lifting Surfaces 

(4) The wing or fin layout routines in the present version of program DEMON3 
can not account for actual gaps between the fin root chords and the body. 
These gaps may be streamwise for undeflected control fins and/or they may 
be due to the fin root chord separated from the surface of the body for a 
deflected fin. Furthermore, the fin deflection angle is only included in the 
flow tangency condition applied at points in the geometrically undeflected 
fin plane. This simplification may cause errors at high fin deflection angles. 
The latest version of the NWCDM-NSTRN program (derived from the code 
described in Reference 8 and extended in connection with on-going work in 
fin tailoring to control the axial location of center of pressure (Reference 29) 
contains subroutines to rotate the control fin in accordance with its 
deflection angle. As such, the influence of fin gap is included in the fin 
loadings as far as in viscid effects are concerned within the limitations of 
supersonic linear theory. Unpublished preliminary results show the same 
trends with gap size and deflection angle as shown by available 
experimental data. The fin tailor version of NWCDM-NSTRN can not 
account for the effects of shocks and shock-boundary layer interactions 
present in the gap. This area of research needs to be studied in detail (with 
Navier-Stokes solvers and suitable grids) since control fin effectiveness and, 
thus, sizing can depend greatly on fin gaps. 


Shock Expansion Method 

(5) As described in an earlier section, program DEMON3 can optionally engage 
shock expansion and Newtonian methods for calculating the pressures 
acting at points on the body meridians and at points along strips on the 
lifting surfaces. The pressures calculated by the method incorporated in 
DEMON3 include effects of aerodynamic interference by modifying the 
flow deflection angles used in the nonlinear pressure calculation methods. 
In the shock expansion procedure, the strip calculations, prior to modifying 
the flow angles, are based on general shock expansion theory. On the 
lifting surfaces, the general approach appears to be adequate. On the body, 
the general method should be upgraded to second-order shock expansion. 
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COMPUTER PROGRAM DEMONS 


General Description 

The DEMON3 computer program, presented in this section, has application as an 
engineering level prediction method to wing-body-tail configurations at angles of 
attack up to 25 deg and for supersonic Mach numbers up to 2.5 if the standard panel 
method/Bernoulli pressure method is used. If the nonlinear (shock 
expansion/Newtonian with interference effects obtained from the panel solutions) 
pressure calculation method is engaged, the upper bound on the Mach number is 
approximately 6 but the body nose and fin leading edge shocks must be attached. On 
the lifting surfaces, the nonlinear pressure calculations compare well with 
experimental data; however, on the forebody the nonlinear pressure calculation on 
the forebody as implemented in the present version of DEMON3 needs to be updated 
to second order. 

The DEMON3 program consists of an executive routine, also designated 
DEMON3, and 79 subroutines. An alphabetical listing of all the routines is provided 
in Table I. The subroutine calling sequence is indicated in Table II. Program 
DEMON3 is written in standard FORTRAN language. Operation requires two (2) 
auxiliary files in addition to the standard input and output units. Core storage for 
execution is approximately 170 kilobytes. On a SunSPARC work station, execution 
time can vary from one minute to one hour depending on the configuration, the 
number of panels selected to model the body and fin sections, and the length of body 
along which body vortices are shed and tracked. 

Program stops are listed in Table HI. Some of the stops print the reason(s) for the 
program stop. All stops specify the subroutine from which they are printed. 

Stepwise Procedure 

A complete configuration consisting of a body with a forward fin set and a tail 
fin set is handled by DEMON3 in four (4) steps. The stepwise procedure is shown in 
Figure 16. The calculations performed for each step are itemized below. 

Step 1. Obtain triplet panel solution for complete body; perform body vortex 
shedding (when applicable) and pressure calculations on the forbody; save 
forebody force and moment data and forebody vortex strength and position 
data for transfer to routines for Step 2. 

Step 2. Obtain constant u-velocity panel solution for the forward fin section under 
the influence of forebody vortices (when present); calculate pressures and 
forces and moments acting on the lifting surfaces and interference shell; 
determine discrete vortex model of the wake for each lifting surface at the 
trailing edge of the forward fin section; transfer wake vortex strengths and 
positions to the routines for Step 3. 
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Step 3. Continue the body shedding calculations along the portion of body 
(afterbody) between the forward fin set and the tail fin set under the 
influence of the vorticity associated with the forward fin set and forebody 
vorticity (when applicable); calculate pressures and forces and moments on 
the afterbody; save vortex strengths and positions and forces and moments 
for inclusion in the final Step4. 

Step 4. Obtain solution for the constant u-velocity panels on the lifting surfaces and 
the interference shell in the tail fin set under the influence of upstream 
vorticity associated with the forebody, forward fin set, and afterbody. 


Input Description 

The following is a description of the input variables required by the present version of 
program DEMONS. The description is given in the sequence of the stepwise procedure 
described above. 


INPUT VARIABLES FOR PROGRAM DEMON3 


Program 

Variable Format Comments 


Input Variables Required for All Rims 

Normally these items are input to Step 1. However, if restarting the 
calculation on the afterbody section with NSTART = 3 in Item 1 below, 
these items are required as input to Step 3. 


Item 1 

(215) 

Beginning and ending steps in calculation. 

NSTART 


Step at which calculations are started or restarted. 
NSTART = lor 3. 

NSTOP 


Step after which calculations are stopped, 
1 < NSTOP < 4. 

Item 2 

(1615) 


NCIR 


Cross Section shape index. 

=0, circular cross section, do not use; use NCIR = 1 
instead 

=1 , elliptic cross section 
=2, arbitrary cross section 
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Program 

Variable 

NCF 

ISYM 

NBLSEP 


NSEPR 

NSMOTH 


NDFUS 

NDPHI 


Format Comments 


Numerical transformation index. 

=0, mapping coefficients are calculated based on body 
shape input and are saved on Unit 8 

=1, mapping coefficients from previous run are read in 
from Unit 8 

Symmetry index. 

=0, right-left flow symmetry 

=1, no symmetry (required if 4> * 0 0 ) 

Body vortex separation index. 

=0, no separation (required if a c = 0 0 and preferred if 
a c > 10 in Item 6) 

=1, laminar separation 

=2, turbulent separation 

Reverse flow separation index. 

=0, no separation, recommended 

=1, laminar separation in reverse flow region 

Vortex-induced velocity smoothing index. 

=0, no smoothing 

=1, smoothing of vortex-induced velocities used in 
pressure calculation 

=2, smoothing of vortex-induced velocities used in 
pressure calculation and smoothing of 8 0/dx 
calculation 

Vortex core model index. 

=0, potential vortex 

=1, diffusion core model; preferred (see RCORE in 
Item 6) 

Unsteady pressure term index. 
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Program 

Variable 


Format 


Comments 


INP 


NXFV 

NFV 

NVP 

NVR 

NVM 

NVA 

NASYM 

Item 3 


=0, omit d<p/dx from C p calculation 
=1, include d<p/dx term; preferred 
Nose force index. 

=0, region between nose tip and XI divided into 10 
evenly spaced sections, and pressure integration is 
performed over this region. Pressures at each 
section are obtained by interpolating between XI 
and the nose tip, which for Bernoulli pressures, is 
assumed to be a stagnation point. Pressures are 
constant over the region and equal to the pressure 
at XI when calculating 2-D nonlinear pressures. 
Preferred. 

=1, no forces calculated ahead of XL Forces at XI are 
input in Item 18. Required value for restart option. 

Number of x-stations at which flowfield is calculated or 
special output generated. See Item 15. 0 < NXFV s 8 

Number of field points for flowfield calculation. See 
Item 16. 0 s NFV < 100 

Number of +T vortices on +y side of body to be input 
for restart calculation. See Item 18. 0 s NVP s 70 

Number of -F reverse flow vortices on +y side of body 
to be input for restart calculation. See Item 20. 
0 s NVR < 30 

Number of -F vortices on -y side of body to be input for 
restart calculation. See Item 21. NVM = 0 if ISYM = 0. 
0 s NVM < 70 


Number of +F reverse flow vortices on -y side of body 
to be input for restart calculation. See Item 22. NVA = 0 
is ISYM = 0. 0 < NVA s 30 

Asymmetric vortex shedding index. See Item 8. 

=0, no forced asymmetry; preferred 

=1, forced asymmetry 

(1215) 
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Program 

Variable 


Format 


Comments 


NHEAD 

NPRNTP 


NPRNTS 


NPRNTV 


NPLOTV 


NPLOTA 


NPRTVL 


Number of title cards in Item 3. NHEAD > 1 
Pressure distribution print index. 

=0, no pressure output except at special x-stations 
specified by Item 15 

=1, pressure distribution output at each x-station 
Vortex separation print index. 

=0, no output 

=1, output at each x-station; preferred 

=2, detailed separation calculation output; for 
debugging purposes only 

Vortex cloud summary output index. 

=0, no vortex cloud output 

=1, vortex cloud output; preferred 

Vortex cloud printer-plot option. 

=0, no plot 

=1, plot full cross-section on a constant scale 
=2, plot upper half cross-section on a constant scale. 

=3, plot full cross-section on a variable scale; preferred. 
Plot frequency index. 

=0, no plots 

=1, plot vortex cloud at x-stations specified by Item 15 

=2, plot vortex cloud at each x-station 

Velocity calculation auxiliary output for debugging 
purposes only (Sub. VLOCTY). 

=0, no output; preferred 

=1, print velocity components at field points; see Items 
15 and 16 
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Program 

Variable 


Format 


Comments 


NCNTR 


NCORE 


NVTRNS 


NRSTRT 


N2DPRB 


Item 4 
TITLE 


=2, print velocity components at body control points 
during pressure calculations. This option can 
produce massive quantities of output (not 
recommended). 

Center vortex control index. 

=0, no center vortex 

=1, center vortex calculated 

Vortex core size index. 

=0, ratio of the local vortex core radius to the local 
body diameter = 0.025, default value 

=1, vortex core, RCORE, input in Item 6. 


Vortex cloud representation index. 

=0, forebody vortices represented by centroids. 

=1, entire set of forebody-shed vortices transferred to 
fin subroutine SPNLD, and thus present during 
afterbody and tail fin calculations. 

Restart control index. 

=0, calculation is not a restart 

=1, calculation is a restart. Additional input required. 
Items 16 and 18. 

Index which controls type of loading calculation 

performed on the body surface. 

=0, Bernoulli pressure coefficients 

=1, shock expansion and Newtonian pressure 
coefficients 

=2, shock expansion and Newtonian pressure 
coefficients corrected with linear theory 

(20 A4) Hollerith information identifying the run 

NHEAD lines of identification. Information to be 

printed at top of output sheets. 
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Program 

Variable 

Format 

Comments 

Item 5 

(8F10.0) 

Reference information used in forming aerodynamics 
coefficients. 

REFS 


Reference area. REFS > 0 

REFL 


Reference length. REFL > 0 

XM 


Xg-coordinate of moment center. 

ZM 


Zg-coordinate of moment center. 

SL 


Body length (L). 

SD 


Body maximum diameter, or equivalent noncircular 
body diameter (D). Used in definition of RE. 

Item 6 

(8F10.0) 


ALPHAC 


Angle of incidence, degrees. 0 0 < a c < 90 0 if a c = 0 0 , 
NBLSEP = 0 in Item 2. 

PHI 


Angle of roll, degrees (0). 

RE 


Reynolds number (VJD/v). 

RCORE 


Ratio of the local vortex core radius to the local body 
diameter. Default value is .025, and maximum 
allowable value permitted in the code is .05*SD/ r 0 . 

XMACH 


Mach number, M„. 

Item 7 

(8F10.0) 



XI 


Initial x-station. XI > 0. 


XF Final x-station. XF > XL 

DX Increment in x for vortex shedding calculation. Typical 

value, DX = D/2 when FDX = 0.0. For variable DX, set 
DX = 0.0 and FDX > 0.0. When N2DPRB < 0 (Item 3), 
DX < XI. 

EMKF Minimum distance of shed vortex starting position from 

body surface. 

=1.0, vortices positioned such that separation point is 
a stagnation point in the crossflow plane 
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Program 

Variable 


Format 


Comments 


RGAM 


VRF 


FDX 


Item 8 (8F10.0) 


E5 


XTABL 


XASYMI 

XASYMF 

DBETA 


>1 .0, minimum radii away from body surface for shed 
vortices. Typical value, EMKF = 1.05. 

Vortex combination factor. 

=0.0, vortices not combined; preferred (set RGAM = 0) 

>0.0 radial distance within which vortices are 
combined. Typical value, RGAM = 0.05 D 

Vortex reduction factor to account for observed decrease 
in vortex strength. 

=0.6, for subsonic flow 

=1 .0, for closed bodies, or supersonic flow; preferred 

Factor for variable x-increments. When FDX > 0.0 
DX = FDX * (local body radius). DX should be input as 
zero when this option is used. 

This item contains only one variable that is of general 
use in program DEMON3 and that is the integration 
error tolerance, E5. The variable XTABL is for 
diagnostic purposes only. The next three variables 
concern the use of forced asymmetry for bodies at very 
high angles of attack. The variables are defined as 
follows: 

Error tolerance for vortex trajectory calculation. 
(Relative error in vortex position.) Typical range, 
E5 = 0.01 to 0.05. Set E5 = 0.0. 

x-location after which a table of corresponding points 
between the real body and the transformed circle is not 
printed. For diagnostic purposes on noncircular shapes. 
Typical value, XTABL = 0.0. 

Initial x-location at which forced asymmetry of 
separation points is used. Typical value, XASYMI = 0.0. 

Final x-location at which forced asymmetry of 
separation points is used. Typical value, XASYMF = 0.0. 

Amount of forced asymmetry for separation points on 
body, degrees. Typical value, DBETA = 0.0. 


Items 9, 10 provide a table of geometry characteristics that must be input for all bodies, 
circular or noncircular in cross-section. 
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Program 


Variable 

Format 

Comments 

Item 9 

(15) 



NXR 


Number of entries in body geometry table 
(1 < NXR < 50). 

Item 10 

(5F10.0) 

NXR lines of input 

XR 


x-stations for geometry table (NXR values, 8 per record). 

R 


r Q , circular body radius at x-stations, or radius of 
transformed circle (NXR values, 8 per record). For an 
ellipse, r Q = (a + b)/2. 

DR 


dr Q /dx, body slope of transformed shape at x-stations 
(NXR values, 8 per record) 

AE 


a, horizontal half-axis of elliptic cross-section (NXR 
values, 8 per record). 

BE 


b, vertical half-axis of elliptic cross-section (NXR values, 
8 per record). 

Items 11 through 14 are included only if NCIR = 2; that is, for a body with arbitrary 
cross-section that must be handled with the numerical transformation and conformal 
mapping procedures. 

Item 11 

(215) 


MNFC 


Number of coefficients describing transformation of 
arbitrary body to a circle (1 < MNFC < 25). 

MXFC 


Number of x-stations at which transformation 
coefficients are defined (1 < MXFC < 10). For a similar 
shape body at all cross-sections, MXFC = 1. 

Item 12 

(8F10.0) 


XFC 


x-stations at which transformation coefficients are 
defined. For a similar cross-section body, MXFC = 1, 
XFC(l) < XI (MXFC values, 8 per record). 


Items 13 and 14 are included when NCF = 0. This block of records is repeated MXFC 
times, once for each x-station input in Item 12. 
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Program 

Variable Format 

Item 13 (15) 

NR 

Item 14 (2F10.0) 

XRC(IJ), YRC(I,J) 


Comments 


Number of coordinate pairs defining the body cross- 
section at the axial station defined by XFC(J). 
(2 < NR < 30). 


Coordinates of right-hand side of body. The convention 
for ordering the coordinates from bottom to top in a 
counter-clockwise fashion, as shown in Figure 17, is 
observed. Right/ left body symmetry is required in the 
numerical mapping. (NR records with one set of 
coordinates per record). 


Items 15 and 16 are included only if NXFV > 0. Item 15 specifies the axial stations at 
which additional output or plots are requested. Item 16 is included only if NFV > 0. Each 
record of this item contains the nondimensionalized y,z-coordinates of a field point at which 
the velocity field is calculated, at each axial station specified in Item 15. 


Item 15 

(8F10.0) 


XFV 


x-station at which field point velocities are calculated or 
at which additional output is printed (NXFV values). 

Item 16 

(2F10.0) 

Omit if NFV = 0 in Item 2. 

YFV, ZFV 


y/r D , z/r 0 -coordinates of field points at which velocity 
field is calculated, expressed as a fraction of the local 
body radius (NFV records with one set of coordinates 
per record). It is important that the field points lie 
outside the body surface. 


Items 17 and 18 are included for a restart calculation only. 

Item 17 (8F10.0) Omit if NRSTRT = 0 in Item 3. 

CPBO(I) Bernoulli pressures at x = XI (input Item 7), for 73 

meridional points from ft = 0 to /3 = 360°; 73 values, 8 
values per record. These pressures should be obtained 
from a previous run containing this x-station. If 
restarting at the beginning of the afterbody section, 
pressures should be taken from values at control points 
on body interference shell constant u- velocity panels. 


Item 18 (6F10.0) Omit if NRSTRT = 0 in Item 3. Item 18 is a single record 

containing the force and moment coefficients at the 
restart point, x = XI. 


- 32 - 



Program 

Variable 

CA 

CN 

CY 

CM 

CR 

CSL 

Item 19 

GAMP 

YP, ZP 

XSHEDP 

Item 20 

GAMR 
XR,ZR 
XSHEDR 
Item 21 


GAMM 


Format 


(4F10.0) 


(4F10.0) 


(4F10.0) 


Comments 

Axial force coefficient. 

Normal force coefficient. 

Side force coefficient. 

Pitching moment coefficient. 

Yawing moment coefficient. 

Rolling moment coefficient. 

Omit if NVP = 0 in Item 2. This Item is a set of NVP 
records. 

r /V„, positive separation vorticity on right side of 
body. 

coordinates of discrete vortices on right side of body at 
starting point (XI). 

x-location at which individual vortex was shed (may be 
0.0). This variable is used only to identify individual 
vortices and permit the user to follow these vortices 
during the calculation. 

Omit if NVR = 0 in Item 2. This item is a set of NVR 
records. It is a convenient way to add an arbitrary cloud 
of vorticity which is to be maintained separately from 
the other body vorticity in the field. 

T/V^ reverse flow or additional vorticity on right side 
of body. 

coordinates of discrete vortices on right side of body at 
starting point (XI). 

x-location at which individual vortex was shed (may be 

0 . 0 ). 

Omit if NVM = 0 or if ISYM = 0 


body fixed coordinate 
system (xg, yg, Zg) 
(Refer to Figure 19) 


This item consists of NVM records. In the case of a 
symmetric flowfield (ISYM = 0), this block of vorticity is 
automatically defined as the mirror image of the 
positive vorticity input in Item 19. 

r/V^ negative separation vorticity on left side of body. 
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Program 

Variable 


Format 


Comments 


YM, ZM 
XSHEDM 


Item 22 (4F10.0) 


GAMA 

YA,ZA 

XSHEDA 


Coordinates of discrete vortices on left side of body at 
starting point (XD- 

x-location at which individual vortex was shed (may be 

0 . 0 ). 


Omit if NVA = 0 or ISYM = 0 


This item consists of NVA records and is analogous to 
Item 20. This block of vorticity is automatically defined 
as the mirror image of the negative vorticity input in 
Item 20 if symmetry is required by ISYM = 0. 

T/V^ reverse flow or additional vorticity on left side of 
body. 

Coordinates of discrete vortices on left side of body at 
starting point (XI). 

x-location at which individual vortex was shed (may be 

0 . 0 ). 


The above 22 items make up the body vortex shedding portion of the input. Items 23 
through 37 are associated with the body triplet panel layout. 


Item 23 (20A4) 

TITLE 

Item 24 (215) 

IXZSYM 


One card of identification information to be printed 
ahead of the paneling output. 


Panel symmetry option. 

=0, panel symmetric half of body for symmetric flow 
(<t> = 0 ) 

=1, panel complete symmetric body using +y side 
geometry for a symmetric body in asymmetric flow 
(<P = 0) 

=-l, panel complete configuration using geometry of 
both sides 


The purpose of IXZSYM is to define the manner in which the body paneling is carried 
out for different flow conditions. A symmetric body in a symmetric flow (0 = 0°) can 
benefit from the IXZSYM = 0 option, which panels only one half of the body and then 
utilizes the symmetry characteristics to reduce computation time. The same symmetric 
body in an asymmetric flow condition ( 0 * 0 °) must be paneled in total, IXZSYM = 1, to 
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Program 

Variable Format Comments 

pick up the nonsymmetry in the loading around the body. The body symmetry is used to 
reduce the input required as described later. 

ITBSYM =0, configuration has top/bottom symmetry 

=1, no top/bottom symmetry (e.g., see Fig. 10) 

Item 25 (413) 

J2 Body geometry specification option. 

= 1 , geometry for arbitrary shaped body 

=-1, circular body defined by cross-section area at XFUS 
stations (not recommended) 

=-2, circular body defined by radius at XFUS stations 

=-3, elliptic body defined by both semi-axes at XFUS 
stations 

J6 =0, cambered body - not available 

=1, body symmetrical with respect to XY-plane 
(uncambered body) 

=-1, circular or elliptic body (preferred value) 

NRADX Number of points used to represent the body segment 

about the circumference. If the configuration is 
symmetric (IXYZSYM = 0, 1), NRADX is input for the 
half section only. If the entire configuration is input 
(IXZSYM = -1), NRADX is input for the full section. The 
half section or full section is divided into NRADX - 1 
equal angles, and the y- and z- coordinates of the panel 
corner points are defined by the intersections of these 
meridian angles with the entire body surface 
(3 < NRADX < 20). 

NFORX Number of axial body stations used to define geometry 

(2 < NFORX < 30). 

The next block of input specifies the geometry of the body for purposes of laying out 
the panels. This input is somewhat redundant with the input in Items 9 through 14; 
however, it is important that the paneling geometry be separate from the previous geometry 
but consistent with it. In many cases, the paneling geometry may be specified in a more 
coarse grid than that required by the separation portion of the program. 
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Program 

Variable 

Format 

Comments 

Item 26 

(8F10.0) 



XFUS 


x-stations at which body paneling geometry is input 
(NFORX values, 10 values per record). 

Item 27 

(8F10.0) 

Included only if J2 = -1 (Item 25). The circular body 
cross-sectional areas are specified at the XFUS stations 
defined in Item 26. 

FUSARD 


Circular body cross-sectional areas at XFUS stations 
(NFORX values, 10 values per record). 

Item 28 

(8F10.0) 

Included only if J2 = -2 (Item 25). The circular body 
radii are specified at the XFUS stations defined in Item 
26. 

FUSRAD 


R, circular body radii at XFUS stations (NFORX values, 
10 values per record). 

Items 29 and 30 are included onlv if T2 = -3 (Item 25). The elliptic body horizontal and 
vertical semi-axes are specified at the XFUS stations defined in Item 21. 

Item 29 

(8F10.0) 


FUSBY 


Elliptic body horizontal semi-axes at XFUS stations 
(NFORX values). 

Item 30 

(8F10.0) 


FUSAZ 


Elliptic body vertical semi-axes at XFUS stations 
(NFORX values). 

Item 31 

(2F10.0) 

Included only if J2 = 1 (Item 25). The y,z-coordinates of 


points on an arbitrary cross-section body are specified in 
Item 31, and there is one set of Item 31 for each of the 
axial stations defined by XFUS in Item 26. The 
convention for ordering the coordinates from bottom to 
top is observed (e.g.. Fig. 10). If the full cross-section is 
to be specified (IXZSYM = -1, Item 24) the ordering 
continues counter-clockwise from the top of the body 
back to and including the bottom point to close the cross 
section. 

YJ, ZJ y,z-coordinates or arbitrary body at XFUS stations 

(NRADX records, 2 coordinates per record). 

Item 31 is repeated NFORX times, once for each XFUS station. 
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Program 

Variable 


Format 


Comments 


The purpose of the following input is to provide a simple means to adjust the actual 
panel layout keeping the geometry specified above unchanged. 

Item 32 (20A4) 


TITLE2 


Item 33 (515) 


IPRINT 


IPRT(l) 


IPRT(2) 

IPRT(3) 


IPRT(4) 


One card of identification information for paneling 

geometry. 

Print control options used to debug specific subroutines. 

Large amounts of output may be generated. 

Panel output option. 

=0, minimum output 

=1, output panel corner coordinates, final source 
strengths, and velocities at panel control point 
(recommended) 

=2, same as IPRINT = 1 and source strengths from 
previous iteration (not recommended) 

=3, print source strength for each iteration (not 
recommended) 

=4, same as IPRINT = 3 and complete aerodynamic 
influence matrix (not recommended) 

=0, no supplementary print 

>0, further triplet panel geometry information from 
subroutines BODPAN, CONFIG, and NEWRAD. 

not used. 

=0, no supplementary print 

=1, print normal influence coefficients from 
subroutines BICOEF and WCOEF. 

=2, print above plus unit velocities, u,v,w, from 
subroutine BICOEF and WCOEF. 

=0, no supplementary print 

>0, print control indices, field point coordinates, and 
triplet panel-induced velocities from subroutine 
BODVEL. 
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Program 

Variable 


Format 


Comments 


>1, print above plus influencing triplet panel strengths 
in subroutine BODVEL. 

>2, print above plus AN arrays from subroutine 
BODVEL. 

>3, print above plus UB, VB, WB arrays from 
subroutine BODVEL. 


If IPRINT < 0, the panel geometry is included with the output specified by I IPRINT I . 
This option provides a useful means to check the paneling arrangement. 


Item 34 

(215) 

These indices provide a means to change the panel 
layout without changing the geometry input. 

KRADX 


Number of meridian lines used to define panel edges on 
the body segment. 

=0, meridian lines are defined by NRADX in Item 29. 

>0, meridian lines are calculated at KRADX equally 
spaced points (3 < KRADX < 20). 



<0, meridian lines are calculated at specified values of 
PHIK in Item 36 (3 < 1 KRADX 1 <20). 

For symmetric configurations (IXZSYM = 0,1)/ KRADX 
is the number of meridians on the half section. For full 
configurations (IXZSYM = -1), KRADX is the number of 
meridians on the full section including meridians at 0 ° 
and 360°. 

KFORX 


Number of axial stations used to define leading and 
trailing edges of panels on the body. 

=0, the panel edges are defined by NFORX and XFUS 
in Items 29 and 30, respectively. 

>0, the panel edges are defined by KFUSK in Item 36 
(2 < KFORX < 30). 

Item 35 

(8F10.0) 

Include only if KRADX < -3 in Item 35. 

PHIK 


Body meridian angles (degrees) with PHIK = 0 0 at the 
bottom of the body and PHIK = 180° at the top 
(1 KRADX 1 values). 

Item 36 

(8F10.0) 

Omit if KFORX = 0 in Item 35. 
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Program 

Variable 

KFUSK 


Format Comments 


Axial stations of leading and trailing edges of body 
panels (KFORX values). XFUSK defines the axial 
positions of the panel edges. If this item is not included, 
the panels are specified by XFUS in Item 26. The 
purpose of this item is to simplify the changing of the 
panel layout without changing the panel geometry. 
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Program 

Variable Format Comments 


Input Variables Required for Step 2 
See notes at end of Namelist $INPUT for variables marked with *. 


Item 1 

(namelist) 

Namelist $INPUT. 

B2 


Exposed fin semispan of horizontal or upper right or 
lower left fins, dimensional. 

B2V 


Exposed fin semispan of vertical or upper left or lower 
right fins, dimensional, default is 0.0. 

CRP 


Rootchord or horizontal or lower left or upper right fins, 
dimensional. 

CRPV 


Rootchord of vertical or upper left or lower right fins, 
dimensional, default is 0.0. 

DELD* 


Deflection angle of vertical lower or upper left fin. 

DELL* 


Deflection angle of horizontal left or lower left fin. 
Positive: trailing edge down, degrees, default is 0.0. 

DELR* 


Deflection angle of horizontal right or upper right fin. 
Positive: trailing edge down, degrees, default is 0.0. 

DELU* 


Deflection angle of vertical upper or lower right fin. 
Positive: trailing edge to right or down, degrees, default 
is 0.0. 

FAC 


FAC = 0.95. Fraction of the constant u-velocity panel 
chord (which contains the centroid) where the control 
point is located. 

FKLE 


Fraction of leading-edge suction converted to normal 
force, default is 0.5. 

FKSE 


Fraction of side-edge suction converted to normal force, 
default is 1.0. 


LVSWP LVSWP = 0 No breaks in fin leading or trailing edges, or 

equal spanwise spacings of panel side 
edges, default value. 

LVSWP h 0 Up to 19 breaks in fin leading or trailing 
edges or up to 19 unequal spanwise 
spacings. 
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Program 

Variable 


Format 


Comments 


MSWD* 

MSWL* 

MSWR* 

MSWU* 

NCW* 

NCWT 

NOUT(I) 


NPR 

NTDAT 

NTPR 

N2DPRF 


Number of spanwise constant u-velocity panels on the 
vertical lower or upper left fin; 1 < MSWD <19, default 
isO. 

Number of spanwise constant u-velocity panels on the 
horizontal left or lower left fin; 1 < MSWL < 19, default 
isO. 

Number of spanwise constant u-velocity panels on 
horizontal right or upper right fin; 1 < MSWR < 19, no 
default. 

Number of spanwise constant u-velocity panels on 
vertical upper or lower right fin; 1 < MSWU < 19, 
default is 0. 

Number of chordwise constant u-velocity panels on the 
fins. 

Number of fin thickness (source) panels in a chordwise 
row, default is 0. Usually NCWT > NCW to capture 
streamwise thickness details on a fin. 

Array of 5 debug print control. Default is NOUT(I) = 0; 
1 < I < 5. These variables produce print in the following 
routines when set equal to 1. 

NOUT(l) Debug print in FINSEC 
NOUT(2) Debug print in RHSFIN 
NOUT(3) Print fin-load results for linear pressures 
NOUT(4) Debug print in BDYSHL 
NOUT(5) Debug print in FINVEL, SPECPR, SPECLD, 
and SPNLD 

Print control option. Large amounts of output may be 
generated if NOUT > 0. Default is 0. 

NTDAT = 0 No thickness input data, default value. 
NTDAT = 1 Fin thickness information to be supplied 
(Item 6) for one or more fins of a finned 
section. 

NTPR = 1 Print debug output from subroutine 
THKVEL, default is 0. 

Index governing type of loading calculation performed 
on fin surfaces. Note: N2DPRF = 1 or 2 not applicable 
to cases including nonzero angular rates. 


- 41 - 



Program 

Variable 


Format 


Comments 


N2DPRF = 0 Linear and Bernoulli pressure coefficients, 
default value. 

N2DPEF = 1 Shock expansion and impact (Newtonian) 
pressure coefficient 

N2DPRF = 2 Shock expansion and impact (Newtonian) 
pressure coefficients corrected with linear 
theory. 

For N2DPEF > 0, further input required (Item 7). 

PHIDIH PHIDIH = 0.0 


PHIFR Dihedral angle of right upper fin, in degrees. Default is 

0.0 for cruciform or planar fins. See sketch below. 

PHIFL Dihedral angle of left lower fin, in degrees, default is 0.0 

for cruciform or planar fins. 

PHIFU Dihedral angle of right lower fin, in degrees. Default is 

90.0 for cruciform fins. 

PHIFD Dihedral angle of left upper fin, in degrees. Default is 

90.0 for cruciform fins. 


Arbitrary Fin Layout 


Cruciform Fin Layout 


Looking Forward 


Zb z b 
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Program 

Variable 


Format 


Comments 


SWLEP* ** 

SWTEP* 

SWTEV* 

THETIT 

THETR~ 

THETL” 

thetlT 

thetd” 

TOLFAC 

XLEBIP 

YBOD(J)* 

ZBOD(J)* 


Vertical or upper left or lower right fin leading-edge 
sweep angle measured in fin planform, positive for 
sweep back, degrees, default is 0.0. 

Horizontal or lower left or upper right fin trailing-edge 
sweep angle measured in fin planform, positive for 
sweep back, degrees, default is 0.0. 

Vertical or upper left or lower right fin trailing-edge 
sweep angle measured in fin planform, positive for 
sweep back, degrees, default is 0.0. 

Polar angle associated with interdigitated or low profile 
four fin layouts, default is 0.0 for cruciform or triform 
fin layouts, 0 < 8 < 90 ° . 

Polar angles of upper right fin, degrees, default is 0.0 for 
cruciform or planar fins. 


Polar angle of lower left fin, degrees, default is 0.0 for 
cruciform or planar fins. 

Polar angle of lower right fin, degrees, default is 90.0 for 
cruciform fins. 


Polar angle of upper left fin, degrees, default is 90.0 for 
cruciform fins. 


Multiplication factor used in the evaluation of the 
tolerance, TLRNC, used throughout the program. 
Default value is 1.0, on VAX 11/750 use 10,000. 

Axial location of fin root chord leading edge measured 
from body nose, dimensional, default is 0.0. 

Zg-coordinate of moment center in body coordinate 
system, feet, default is 0.0. 

yg-coordinate of rootchord leading-edge, 1 < J < 4. User 
must specify a value for each wing or fin present in 
input description. 

Zg-coordinate or rootchord leading edge, 1 < J < 4. User 
must specify a value for each fin present in input 
description. 


* For fins with breaks in leading and/ or trailing-edge sweep, use default values. 

** These angles specified only for arbitrary fin layouts, refer to preceding sketch. 
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The following relations must hold: 


1. NWBP = NPANLS + NCWB*NBDCR < 300. 

NPANLS = NCW * (MSWR + MSWL + MWSU + MSWD) <150. 

NBIP < 150 (NBIP = number of panels in interference shell.)*** 

2. Also, MSWR, MSWL, MSWU, and MSWD should be at least five for valid fin training- 
edge vorticity characteristics. 

3. When running symmetric case, set PHI = 0.0 in Item 6, IXZSYM = 0 in Item 33 for all 
runs (preceding Step 2 input). If fins are deflected symmetrically, set MSWL, MSWU, 
MSWD, and DELL, DELU, DELD as follows: 

(1 ) for cruciform fins, MSWR + 0, DELR + 0.0., 

set MSWL = 0 DELL = DELR 

MSWU = 0 DELU = 0 

MSWD = 0 DELD = 0 

(2) for arbitrary but symmetric fin layout, MSWR * 0, MSWU + 0, 2 cases: 

a no fin deflection 
set MSWL = 0 
MSWD = 0 

b with fin deflection, set DELR * 0.0, DELU * 0.0, 
set MSWL = MSWU DELD = DELR 

MSWD = MSWR DELL = DELU 

The following four variables are used below in the description of the input variables. 
The terms "right" and "left" refer to an observer looking forward. 

MSWRP: Number of chord wise rows on horizontal right or upper right fin + 1; upper 
right, lower left fin, etc., terminology is used for arbitrary fin layouts shown 
in preceding sketch. 

MSWRP = MSWR + 1 

MSWLP: Number of chord wise rows on horizontal left or lower left fin + 1; 

MSWLP = MSWR + 1 

MSWUP: Number of chordwise rows on vertical upper or lower right fin + 1; 

MSWUP = MSWU + 1 

MSWDP: Number of chordwise rows on vertical lower or upper left fin + 1; 

MSWDP = MSWD + 1 


*** NBIP is set by user. It is determined by selection of triplet panel layout on the body over 
the length of the fin rootchord. 
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Program 

Variable 


Format 


Comments 


Item 2 (3F10.0) 

YRT(KJ) 

VSWLER(KJ) 

VSWTEROCJ) 

Item 3 (3F10.0) 

YLTCKJ) 

VSWLEL(KJ) 

Item 4 (3F10.0) 

ZUT(KJ) 

VSWLEUCKJ) 

VSWTEUOCJ) 


Optional input for fin-body combination (if LVSWP = 0). 

Distance from fin rootchord to the constant u- velocity 
panel outboard side edge on right horizontal or upper 
right fin, dimensional 1 < KJ < MSWRP, (MSWRP < 20), 
YRT(l) = 0.0, YRT(MSWRP) = B2. 

Leading-edge sweep of fin between YRT(KJ - 1) and 
YRT(KJ), positive for sweep back, degrees, 
1 < KJ < MSWRP, (MSWRP s 20), VSWLER(l) = 0.0. 
Sweep angle is measured in fin planform plane. 

Trailing-edge sweep of fin between YRT(KJ - 1) and 
YRT(KJ), positive for sweep back, degrees, 
1 s KJ s MSWRP, (MSWRP s 20), NSWTER(l) = 0.0. 
Sweep angle is measured in fin planform plane. 

Optional input for fin-body combination (if LVSWP * 0). 
Omit if MSWL = 0. 

distance (negative) from fin rootchord to the constant u- 
velocity panel outboard side edge on left horizontal or 
lower left fin, dimensional, 1 < KJ s MSWLP, 
(MSWLP < 20), YLT(l) = 0.0, YLT(MSWLP) = -B2. 

Leading-edge sweep of fin between YLT(KJ - 1) and 
YLT(KJ), negative for sweep back, degrees, 
1 s KJ < MSWLP, (MSWLP < 20), VWSTEL(l) = 0.0. 
Sweep angle is measured in fin planform plane. 

Optional input for cruciform fin-body combination (if 
LVSWP * 0). Omit if MSWU = 0. 

Distance from fin rootchord to the constant u- velocity 
panel outboard side edge on upper vertical or lower 
right fin, dimensional, 1 < KJ < MSWUP, (MSWUP < 20), 
ZUT(l) = 0.0, ZUT(MSWUP) = B2V. 


Leading-edge sweep of fin between ZUT(KJ - 1) and 
ZUT(KJ), positive for sweep back, degrees, 
1 < KJ < MSWUP, (MSWUP < 20), NSWLEU(l) = 0.0. 
Sweep angle if measured in fin planform plane. 

Trailing-edge sweep of fin between ZUT(KJ - 1) and 
ZUT(KJ), positive for sweep back, degrees, 
1 < KJ < MSWUP, (MSWUP < 20) VSWTEU(l) = 0.0. 
Sweep angle is measured in fin planform plane. 
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Program 

Variable Format 


Item 5 (3F10.0) 


ZDT(KJ) 


VSWLED(KJ) 


VSWTEEKKJ) 


Item 6 


Item 6(a) (1015) 


MSWT 


Comments 

Optional input for cruciform wing-body combination (if 
LVSWP x 0). Omit if MSWD = 0. 

Distance (negative) from fin rootchord to the constant u- 
velocity panel outboard side edge on lower or upper left 
fin, dimensional, 1 < KJ < MSWDP, (MSWDP < 20), 
ZDT(1 ) = 0.0, ZDT(MSWDP) = -B2V. 

Leading-edge sweep of fin between ZDT(KJ - 1) and 
ZDT(KJ), negative for sweep back, degrees, 
1 < KJ < MSWDPL (MSWDP < 20), VSWLED(l) = 0.0. 
Sweep angle is measured in fin planform plane. 

Trailing-edge sweep of fin between ZDT(KJ - 1) and 
ZDT(KJ), negative for sweep back, degrees, 
1 s KJ < MSWDP, (MSWDP s 20), VSWTED(l) = 0.0. 
Sweep angle is measured in fin planform plane. 

Optional thickness input data when NTDAT * 0 in 
namelist $INPUT. The maximum number of source 
panels for one finned section cannot exceed 400; 
therefore, 

f ? S (MSWT * NCWT) s 400; 

1 

NCWT is specified in namelist $INPUT. 

Information in Items 6(a), 6(b), and 6(c) is read in by 
subroutine THKIN for the right horizontal fin. 

Number of source panels in the spanwise direction, 
1 < MSWT s 19. 


LVSWT LVSWT = 0 No breaks in fin leading or trailing edges, or 

equal spanwise spacings of source panel 
sides, default is 0. 

LVSWT = 1 Up to 29 breaks in fin leading or trailing 
edges or up to 19 unequal spanwise 
spacings. 

NUNIS NUNIS = 0 Chordwise thickness distribution varies 

over the span. 

NUNIS = 1 Chordwise thickness distribution constant 
over the span. 

Item 6(b) (3F10.0) Optional input for LVSWT = 1. 
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Program 

Variable 


Format 


Comments 


YTH(IJ) 

SWLET(IJ) 

SWTET(IJ) 

Item 6(c) 

THET(K) 

Item 6(d) 
Item 6(e) 
Item 6(f) 

Item 7 

Item 7(a) 
CONSTK 


Distance from fin chord to the source panel outboard 
side edge, 1 < J < MSWT + 1. 

Leading-edge sweep of fin between YTH (1,J-1) and 
YTH(IJ), positive for sweep back, degrees, 
1 < J s MSWT + 1. 

Trailing-edge sweep of fin between YTH(1,J-1) and 
YTH(1,J) positive for sweep back, degrees, 
1 < J < MSWT + 1. 

(8F10.0) Optional input specifying streamwise thickness slopes 

read in by subroutine THETIN in groups of NCWT 
values. 


NUNIS = 1 K = 1, NCWT. Input one set of NCWT 
values. 

NUNIS = 0 K = 1, (NCWT*MSWT). Input MSWT sets, 
each set containing NCWT values. 
1 < K < 400. 

(8F10.0) Optional thickness input for left fin for cases with 

nonzero roll or asymmetric fin deflection (MSWL > 0). 
All input same as for right fin above. Items 6(a), 6(b), 
and 6(c). Refer to Item 3 above for positive directions. 

(8F10.0) Optional thickness input for upper fin when MSWU > 0 

in namelist $INPUT. Same input as for right fin. Items 
6(a), 6(b), and 6(c). Refer to Item 4 above for positive 
directions. 

(8F10.0) Optional thickness input for lower fin when MSWD > 0 

in namelist $INPUT. All inputs same as for right fin. 
Items 6(a), 6(b), and 6(c). Refer to Item 5 above for 
positive directions. 

Optional input for calculation of 2-D nonlinear 
pressures on fins. Read when N2DPRF > 0. 


(F10.0) 


Constant used in Newtonian pressure coefficient, 
normally CONSTK = 2. 
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Program 

Variable 

Item 7(b) 
NSEG 

Item 7(c) 
THETA(J) 

Item 7(c) 
THETB(J) 


Format Comments 

(15) 


Number of 2-D segments used to describe fin profile, 
2 s NSEG s 100. 


(8F10.0) 


Slope angles of 2-D segments on upper surface of fin 
measured relative to fin chordal plane, degrees, 
1 < J < NSEG. 


(8F10.0) 


Slope angles of 2-D segments on lower surface of fin 
measured relative to fin chordal plane, degrees. Note: 
Items 7(b) through 7(d) are repeated for each chordwise 
row of control points on each fin. That is, there are 
MSWR + MSWL + MSWD + MSWU repetitions of Items 
7(b) through 7(d). 
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Input Variables Required for Step 3 

Optional input if NSTART = 1 and NSTOP > 2 in Item 1 for all runs. The following 
items from the input required for all runs are repeated here as input to Step 3 for the 
afterbody calculation. 

Item 1 Defined as Item 2 in Step 1 input section. 

Item 2 Defined as Item 3 in Step 1 input section. 

Item 3 Defined as Item 4 in Step 1 input section. 

Item 4 Defined as Item 7 in Step 1 input section. 

Input Variables Required for Step 4 

Optional input if NSTOP = 4 in Item 1 for all runs. Same as for Step 2, applied to the 
tail section. 


Sample Cases, Input Description 

This section describes a sample case input file corresponding to the model of 
Reference 23 and shown in Figure 8. The configuration has a monoplane wing and tail fin 
sections and a noncircular body. The body of the configuration is nearly circular over the 
ogive nose then transitions to an elongated (vertical) circular section on the forebody. 
Before the wing leading edge, the cross section transitions to have a circular bottom and a 
rectangular top. Over the wing root chord, the cross section transitions back to an elongated 
(vertical) circle and maintains this shape to the base. Tail fins are located at the rear of the 
configuration. 

The input file for sample case 1 is shown in Figure 18. Input for calculation step 1 is 
shown on pages 1 through 10, Figure 18(a) through 18(j). Input for calculation steps 2 
through 4 follow on pages 11 through 13, respectively. The input will now be described in 
detail. 

NSTART = 1 and NSTOP = 4 in input item 1 indicate that the calculation procedure will 
start at step 1 (forebody) and stop after step 4 (second fin section). 

Input item 2 consists of many control variables for the body calculation procedure. In 
input item 2, NCIR = 3 indicates that the body cross section is noncircular and changes 
shape (contour) along the length of body. NCF = 1 indicates that the conformal mapping 
coefficient for the noncircular body cross sections have been calculated previously and are 
read from a data file. ISYM = 0 indicates that the geometry and flow conditions have right- 
left symmetry (symmetric geometry and 0 = 0°). NBLSEP = 1 indicates that the Stratford 
laminar separation criteria is used to determine if body separation occurs. No reverse flow 
separation is considered, NSEPR = 0, and vortex induced velocity smoothing is not done, 
NSMOTH = 0. NDFUS = 1 and NDPHI = 1 indicate that a vortex core model is used and 
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that the dQ/dx term is included in the pressure calculation. INP = 0 indicates that nose 
forces are estimated ahead of the first axial station XL The remaining variables in input 
item 2 are zero. These variables govern special output stations, field point velocity 
calculations, and input vortices for a restart calculation. These options are not invoked for 
this sample case. 

Input item 3 consists of additional control variables for the body calculation 
procedure. NHEAD = 1 indicates that one. title line for the run is to be read in item 4. 
NPRNTP, NPRNTS, and NPRNTV are all set to 1 indicating that the pressure distribution, 
separation summary, and vortex summary output pages will be printed at each axial station 
along the body. NPLOTV = 3 indicates that the vortex cloud plots will be plotted on a 
variable scale from one axial station to the next. NPLOTA = 2 indicates that a vortex cloud 
plot will be made at each axial station at which vortices are present. NPRTVL = 0 indicates 
that no debug output is requested from the velocity calculation subroutine VLOCTY. 
NCNTR = 0 indicates that the center vortex is not included in the calculation. NCORE = 0 
indicates that the default vortex core size of 0.025 of the local body diameter is used. 
NVTRNS = 1 indicates that the entire vortex field is transferred to the fin section calculation 
procedure of step 2. NRSTRT = 0 indicates that the calculation is not a restart. N2DPRB = 0 
indicates that Bernoulli pressure coefficients are calculated. 

Input item 4 is a single line (NHEAD) describing the run. 

Input item 5 specifies the reference area, REFS = 159.9984, the reference length, 
REFL = 9.176, the axial coordinate of the moment center, XM = 25.15, the vertical location of 
the moment center, ZM = 0.0, the body length =38.462, and the maximum body diameter, 
SD = 3.0. 

Input item 6 contains the flow conditions for the run. The angle of attack is 10.0° , 
ALPHAC, the body roll angle is 0.0°, PHI, the Reynolds number based on diameter 
is 1.26xl0 8 , RE, the vortex core size is 0.025, RCORE, and the freestream Mach number is 1.6, 
XMACH. 

Input item 7 contains variables associated with the calculation procedure along the 
forebody. The first axial station, XI, is 1.0, the final station, XF, is 18.784, and the axial 
marching step size, DX, is 1.0. The final axial station corresponds to the wing root chord 
leading edge location. EMKF = 1.05 indicates that vortices are placed a minimum of 1.05 
times the local radius away from the body. RGAM = 0 indicates that vortices are not 
combined. The vortex reduction factor, VRF, is set to 1.0. FDX = 0.0 indicates that a variable 
axial step size is not used. 

In input item 8 the tolerance for the vortex tracking integration subroutine is set to 0.05. 

NXR = 14 in input item 9 indicates that 14 points make up the body radius table input 
in item 10. Item 10 is a set of NXR input lines specifying the body radius distribution table. 
For noncircular bodies, the radius corresponds to transformed radius from the conformal 
mapping procedure. The axial station, XR, the transformed radius, R, and the "body slope" 
d(R) / dx are input. 
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MNFC = 25 and MXFC = 6 in input item 11 indicate that 25 mapping coefficients are 
used to model the noncircular cross sections at 6 axial stations. The program interpolates at 
axial stations between the 6 input stations. Item 12 specifies the 6 axial stations where the 
body cross section is input. Items 13 and 14 specify the cross section shape and are included 
for each of the 6 stations specified in item 12. Item 13 specifies the number of coordinate 
points specifying the cross section, NR. Item 14 specifies the coordinates, XRC(I) and 
YRCCD, of the right half of the cross section. Right/left symmetry is required. 

Input items 15 through 22 are not required by this sample case. Most of these items 
pertain to a restart calculation. Items 23 through 36 pertain to the body paneling. 

Input item 23 is a single title line for the paneling part of the input. IXZSYM = 0 in 
item 24 indicates that the body and flow are symmetric and that only half the body will be 
paneled. ITBSYM = 1 indicates that the body does not have top/bottom symmetry. In 
item 25, J2 = 1 indicates the body cross sections are noncircular and nonelliptic, and J6 = -1, 
NRADX = 15, and NFORX = 24 indicate that 14 panels around the half circumference by 23 
panels along the body length are used to model the body. This results in a total of 322 
panels on the half body. Only half of the body is modeled because of symmetry. 

Item 26 specifies the 24 (NFORX) axial stations specifying the panel edge points, XFUS. 
The first point is at 0.0 and the last point is at the end of the body 38.462. Items 27 
through 30 are not included because the model is not circular or elliptic. Item 31 contains 
the coordinates of the panel edge point around the circumference of the body cross sections. 
Item 31 is repeated for each of the 24 XFUS axial stations. Item 31 contains the 15 (NRADX) 
coordinates YJ and ZJ for each axial station XFUS. 

Input item 32 is a single input title line for the additional paneling options in the 
program. In item 33, IPRINT = 1 indicates that panel comer points, final source strengths, 
and velocities at panel control points are printed. The following variable in item 33, IPRT(l) 
through EPRT(4), are set to 0 indicating that no additional supplementary (or debug) output 
is requested. KRADX = 0 and KFORX = 0 in item 24 indicates that the paneling geometry 
defined in items 23 through 31 are to be used to panel the body. Items 35 and 36 are not 
required for this sample case. 

This completes the input description for the required variables for step 1 of the 
calculation procedure. The variables for step 2 through 4 are described next. 

The variables for step 2 of the calculation procedure are input through NAMELIST 
$INPUT shown in Figure 18(k). XLEBIP = 18.784 indicates that the body interference shell 
(also the wing leading edge) start at axial station 18.784. The right wing is attached to the 
body at coordinates YBOD(l) = 1.5 and ZBOD(l) = 0.0, and the left wing is attached at 
coordinates YBOEK2) = -1.5 and ZBOD(2) = 0.0. The exposed span, B2, is 8.5, the root chord, 
CRP is 11.732, and the leading edge sweep, SWLEP, is 44.0 °. A 6x6 chordwise by span wise 
panel layout is used on the wing, NCW x MSWR. Wing thickness is not modeled for this 
case, NTDAT = 0. NOUT(5) = 1 indicates that detailed output will be printed from the fin 
load calculation procedure. 

The variables for step 3 of the calculation procedure are shown in Figure 18(1). The 
values required by step 3 are the same as items 2, 3, 4, and 7 of step 1. Note that the axial 
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starting location, XI, is set to the wing root chord trailing edge location, and that XF is set to 
the tail leading edge location. 

The variables for step 4 of the calculation procedure are shown in Figure 18(m). These 
variables are same as those of step 2; only they specify the geometry and paneling 
information for the tail fins. 

This completes the sample case input description. A brief description of the output for 
this sample case follows. 


Sample Cases, Output Description 

This section describes a sample case output file corresponding to the model of 
Reference 23 and shown in Figure 8. The configuration has a monoplane wing and tail fin 
sections and a noncircular body. The body of the configuration is nearly circular over the 
ogive nose, then transitions to an elongated (vertical) circular section on the forebody. 
Before the wing leading edge, the cross section transitions to have a circular bottom and a 
rectangular top. Over the wing root chord, the cross section transitions back to an elongated 
(vertical) circle and maintains this shape to the base. Tail fins are located at the rear of the 
configuration. The input file for this sample case is shown in Figure 18. 

Figure 19 depicts the coordinate systems and the force and moment sign conventions 
used in DEMON3. The normal force is positive up, the side force is positive to the right, the 
pitching moment is positive nose up, and the yawing moment is positive nose to right. 

The sample case output file is shown in Figure 20. The program title page is output 
page 1, Figure 20(a). Output for step 1 of the calculation procedure (forebody) is shown on 
pages 2 through 95, Figure 20(b) through 20(v). Output for step 2 of the calculation 
procedure (wing section) is shown on pages 96 through 119, Figure 20(w) through 20(aq). 
Output for step 3 of the calculation procedure (afterbody) is shown on pages 120 
through 131, Figure 20 (ar) through 20 (aw). Output for step 4 of the calculation procedure 
(tail fin section) is shown on pages 132 through 153, Figure 20(ax) through 20(bp). A 
summary of the total forces and moments acting on the vehicle are shown on page 154, 
Figure 20(bq). The following is a general description of the sample case output file. Details 
will be provided for calculation steps 1 and 2. Steps 3 and 4 are similar to steps 1 and 2, 
respectively, and will be described briefly. 

Output page 2 displays the flow conditions and prints the program control variables 
for step 1. Page 3 contains the input geometry table, x, r, and dr/dx. Output pages 4 
through 7 contain the numerical conformal mapping coefficients for the six noncircular 
body cross sections. 

Pages 8 through 13 contain the body paneling geometry information. Additional 
paneling information is shown on page 14. The panel comer points are printed on pages 15 
through 20, and the body centroid points are shown on Pages 21 through 26. The panel 
areas and inclination angles are shown on pages 27 through 32. The panels source solution 
for a c = 0 ° is printed on pages 32 and 33. Surface velocities on the body are shown on 
pages 33 through 38. The panel solution at angle of attack is printed on pages 38 and 39. 
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Output pages 40 through 95 contain the pressure distribution, vortex shedding, and 
vortex tracking output for the forebody. Page 40 and 41 contain the pressure distribution 
and sectional loads at the first axial station, x = 1.0. Page 42 is a summary of separation 
information; it is seen that separation does not occur at this axial station. Pages similar to 
pages 40 through 42 are printed as the procedure marches down the body. The pressure 
distribution at the first axial station where separation occurs (x = 7.0) is shown on pages 52 
and 53. The separation summary information is shown on page 54. The marching 
procedure continues down the body to the wing leading edge. A summary of the vortex 
field at the wing leading edge is shown on page 90. The pressure distribution at this station 
is printed on pages 91 and 92. A line printer plot of the vortex field is depicted on page 93. 
New vortices shed at this station are shown on page 94. Page 94 and 95 contain a summary 
of the total loads and the vortex field at the end of the forebody. 

Page 96 begins the output for the wing section force calculation, step 2. The 
NAMELIST $INPUT is printed on this page. Page 97 contains a summary of the wing 
planform geometry, sweep angles, and attachment points to the body. Paneling information 
for the body interference shell is printed on page 98. 

Velocities induced at the wing control points by the body triplet panels are shown on 
page 99. Vortex-induced velocities at the panel control points are shown on page 100. 
Velocities from the constant u-velocity panel modeling the wing and body interference shell 
at wing control points are printed on page 102. 

The pressure distribution at points on the body interference shell are shown on pages 
103 and 104. The results of integration of the pressure distribution are shown on pages 104 
through 106. A summary of the sectional loads on the body interference shell is shown on 
pages 106 and 107. 

The total velocities at the fin control points are shown on page 108, and the velocities at 
panel outboard edges are printed on page 109. Velocities and pressures on the upper and 
lower wing surfaces are shown on page 111. The pressure loadings at these points is shown 
on page 112. A summary of the linear pressure wing loads is printed on page 113. 

Pages 114 and 115 contain the wing spanwise load distribution information used to 
determine the leading edge and side edge suction forces on the wing. The wing vorticity 
from these forces is also determined. 

Page 116 Contains a summary of the Bernoulli pressure loads acting on the wing. 

Pages 117 through 119 contain the spanwise load distributions for the Bernoulli 
loadings. This information is used to determine the vorticity shed from the wing. 

Page 120 begins the output for the afterbody calculation, step 3. This procedure is 
similar to the forebody calculation (step 1) except the procedure starts with initial vorticity 
from the forebody and wings. The pressure distribution on the body is calculated, and 
vorticity is shed and tracked along the afterbody from the wing trailing edge to the tail 
leading edge. The output is shown on pages 120 through 131. These pages are similar to the 
forebody pages described earlier. 
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The output for step 4 of the calculation procedure (tail fin section) is shown on pages 
132 through 153. This output is similar to the output described for step 2 of the calculation 
procedure (wing section). 

The final output page is number 154 which contains a summary of the forces and 
moments acting on the vehicle. Forces and moments are given in both the body rolled and 
unrolled coordinate systems as depicted in Figure 19. 


CONCLUSIONS 


Program DEMON3 for the engineering level prediction of pressures, component loads, 
and overall forces and moments acting on supersonic configurations is described. The 
addressable configuration consists of a body which may have noncircular cross section and 
up to two fin sets. Within a fin set, the lifting surfaces may be in cruciform, trifrom, planar, 
and low profile layouts. The planform of the lifting surfaces may have breaks in sweep, and 
are not limited to any restrictions in aspect ratio or taper ratio. Program DEMON3 is a 
combination of proven methodology including low order panel methods, body vortex 
shedding and tracking schemes, and optional nonlinear pressure calculation methods 
including effects of interference. The particular methods allow application to 
nonconventional configurations at angles of attack up to 25 deg, for Mach numbers up to 
approximately 2.5 if the linear theory/Bemoulli pressure methods are used, or up to Mach 6 
if the optional nonlinear shock expansion and/or Newtonian pressure methods are 
engaged. 

Comparisons with experimental data and other predictions for realistic, 
nonconventional configurations show that the DEMON3 program is capable of predicting 
overall aerodynamic forces, moments, and pressure distributions to a fair to good degree of 
accuracy. However, the nonlinear pressure calculation on the body should be extended to 
second order shock expansion. 

On the basis of the findings reported herein, the following recommendations are 
offered: 

• Additional testing needs to be performed to assess the range of applicability of the 
program further. 

• Program DEMON3 can economically generate aerodynamic force distributions acting 
on the body and the lifting surfaces of nonconventional configurations. This capability 
makes DEMON3 a useful tool in preliminary structural analysis of the components of 
the configuration. In this application, DEMON3 should be equipped with a post 
processor for conversion of forces at aerodynamic control points to forces at user 
supplied structural analysis node or grid points. This technology is presently available 
(Refs. 8 and 12). 
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• The treatment of the center vortex in the vortex tracking scheme should be updated with 
the methodology very recently developed at NEAR for other aerodynamic prediction 
methods. 
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TABLE! 


Routines of program DEMON3 
The designation .f signifies routines in source code 


1. aicbod.f 

2. aicwgb.f 

3. amap.f 

4. asum.f 

5. bdyshl-f 

6. bipan.f 

7. bmap.f 

8. bodpan.f 

9. bodvel.f 

10. cmap.f 

11. combin.f 

12. config.f 

13. confor.f 

14. cpbem.f 

15. demon3.f 

16. dfeqkm.f 


17. dphidti 

18. dtable.f 

19. dzdnu.f 

20. edges, f 

21. elnose.f 

22. f.f 

23. fincon.f 

24. finsec.f 

25. finvel.f 

26. forcep.f 

27. formom.f 

28. fpvelf 

29. gaussj.f 

30. geom.f 

31. inpt.f 

32. layout.f 


33. lubksb.f 

34. ludcmp.f 

35. matmul-f 

36. mxout.f 

37. newrad.f 

38. panel.f 

39. pinpnlf 

40. plotv.f 

41. prcor.f 

42. prcorb.f 

43. prstrt.f 

44. rhsfin.f 

45. rotbw.f 

46. rotfw.f 

47. rotwb.f 

48. rotwf.f 


49. seprat.f 

50. shape.f 

51. shkaglf 

52. shkexp.f 

53. shkxb.f 

54. smooth.f 

55. solve.f 

56. sorpan.f 

57. specld.f 

58. specpr.f 

59. spnld.f 

60. srfvel.f 

61. sum.f 

62. suppan.f 

63. thdih.f 

64. thetax.f 


65. thetin.f 

66. thkin.f 

67. thklyt.f 

68. thkout.f 

69. thkvel.f 

70. tolds.f 

71. vcentr.f 

72. velnor.f 

73. velo.f 

74. veloth.f 

75. vlocty.f 

76. volume.f 

77. vtable.f 

78. xlocatf 

79. z.f 
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TABLE II 


SUBROUTINE CALLING CHART - PROGRAM DEMON3 


PROGRAM 

DEMON3 

l-TOLDS 

I — INPT — PRSTRT 

I I -CONFOR — Z 

I --ELNOSE — F — SHAPE 

I I l-BMAP — AMAP — Z 

I I l—CMAP — DZDNU 

I I I I-Z 

I I I— AMAP — Z 

I I i -BODVEL— SUPPAN 

I I I I — SORPAN 

I I l-VLOCTY— SUM —DZDNU 

I I I -SMOOTH 

I I l-ASUM 

I I -SHAPE 

I I -BMAP —AMAP — Z 

I I — CMAP —DZDNU 

I I l-Z 

I I -GEOM —CONFIG 

I I I— NEWRAD 

I I I -BODPAN— PANEL 

I I — AICBOD — BODVEL — SUPPAN 

I I I -SORPAN 

I I —SOLVE — LUDCMP 

I I I— LUBKSB 

I I I— MATMUL 

I I I -VOLUME 

I I — SRFVEL 

I l-VLOCTY — SUM —DZDNU 

I I I— SMOOTH 

I I l-ASUM 

I I -DPHIDT— SMOOTH 

I I -THETAX — SHAPE 

I I l-BMAP — AMAP — Z 

I I-SHKXB — THDIH 

I I — PRCORB 

I I— FORCEP 

I l-PLOTV — PLOTA2 — PLOTA8 — PLOTA6 

I I I I — PLOTA5 

I I I — PLOTA7 

I I -FPVEL —CMAP —DZDNU 

I I I l-Z 

I I I -BODVEL— SUPPAN 

^ j | I -SORPAN 
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T T T--VLOCTY — SUM — DZDNU 

I I I -SMOOTH 

I I l-ASUM 

I I -SEPRAT 

I I -VIABLE 

I I -DFEQKM — F —SHAPE 

I I l-BMAP — AMAP — Z 

I I l-CMAP — DZDNU 

I I I l-Z 

I I I -AMAP — Z 

I I I -BODVEL— SUPPAN 

I I I l-SORPAN 

I I l-VLOCTY— SUM —DZDNU 

I I I -SMOOTH 

I I l-ASUM 

I I -AMAP — Z 

I I — DTABLE 

i I — VCENTR 

I I — COMBIN 

I l-TOLDS 

I -FINCON— FINSEC — LAYOUT 
I l-THKIN — EDGES 

I 1 I -THETIN 

I I -THKOUT 

I I — THKLYT 

I I — XLOCAT 

I I -BIPAN 

I -AICWGB— VELNOR— ROTWF 
I I — ROTWB 

I l-VELO 

I l-ROTFW 

I l-ROTBW 

I -BODVEL— SUPPAN 
I l-SORPAN 
I — CMAP — DZDNU 
I l-Z 

l-VLOCTY— SUM —DZDNU 
I I -SMOOTH 

I I -ASUM 

I -RHSFIN — THKVEL— PINPNL 
I I I -ROTWF 

I I I— VELOTH 

I I I -ROTFW 

I I -LUDCMP 
I I— LUBKSB 

I l-MATMUL 
I — BDYSHL — CPBERN 
I I -FORMOM 
| |-TOLDS 
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-FINVEL — VELNOR — ROTWF 
I I —ROTWB 
I I -VELO 
I I — ROTFW 

I I — ROTBW 

I -ROTWF 
I -ROTFW 

-SPECPR-— ROTWF 

I -SHKEXP— ROTWF 
I I -SHKAGL 
I -VELNOR — ROTWF 
I I -ROTWB 
I I -VELO 
I I -ROTFW 
I I -ROTBW 
I — THKVEL — PINPNL 
I I -ROTWF 
I I -VELOTH 
I I -ROTFW 
I -PINPNL 
I — CPBERN 
I — PRCOR 

I -SPECLD— ROTWF 
I -ROTFW 
I -TOLDS 

l-SPNLD — ROTWF 
I -ROTFW 
I -TOLDS 


PROGRAM 

END 
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TABLE III 


PROGRAM STOPS 


stop number 

explanation 

subroutine 

20 


bodpan 

50 


bodpan 

60 


bodpan 

unnumbered 

normal stop 

demon3 

201 


dtable 

202 


dtable 

unnumbered 

bad argument 

dtable 

40 


elnose 

41 


elnose 

44 


elnose 

45 


elnose 

46 


elnose 

5 


inpt 

unnumbered 

too many vorticies 

inpt 

unnumbered 

np dimension to high 

ludcmp 

unnumbered 

singular matrix 

ludcmp 

120 


newrad 

130 


newrad 

90 


pinpnl 

unnumbered 

illegal argument for sqrt 

shkagl 

unnumbered 

illegal argument for acos 

shkagl 

unnumbered 

illegal argument for asin 

shkxb 

unnumbered 

hs is negative 

shkxb 

unnumbered 

to many triplet panels 

solve 

201 


vtable 

202 


vtable 

unnumbered 

bad argument 

vtable 

unnumbered 

cannot find triplet panel for interference shell 

xlocat 
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TWO DIMENSIONAL FLOW 


/ 



CONSTANT SOURCE 

(a) 


CONSTANT VORTEX CONSTANT TRIPLET 

( LINEARLY VARYING (UNl-OIRECTIONAl) 

DOUBLET SHEET) 

Supersonic triple! concept. 



(b) Triplet panel grouping on arbitrary body at incidence. 


Figure 1. - Triplet panel concept and layout on 
the body. 
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Figure 2. - Coordinate systems and typical panel layout 
shown for one fin and quarter of the inter- 
ference shell. 
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Figure 4. - Fin wake model. 





Figure 5. 


Span loading and vorticity distribution on a delta 
wing in the presence of an external vortex. 
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the canard-body-tail model in Reference 22, M„, = 2.5, 0 = 0 °, and 
So 4 = ±5° roll control. 
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Figure 9. - Concluded. 
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160.0 irr 
9.1 76 in 
25.15 in 


DEMON3 w/BDYSHD 
DEMON3 w/o BDYSHD 
MISL3 (mod) 

MISDAT (Ref. 28) 
Experiment (Ref. 23) 


(a) Normal-force coefficient 



(b) Pitching-moment coefficient 


Figure 10. - Comparison of predicted and measured aerodynamic characteristics for the 
wing-body-tail model in Reference 23, M = 2.0 and $ = 0*. 
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Paneling layout 


Figure 13. - Panel and pressure tap layout on a rectangular 
wing. 
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Figure 14.— Theoretical and experimental chordwise pressure distribution 
on a rectangular wing at M eo =2.86 and a c =10.3 deg. 
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Figure 15.— Pressure distributions acting on upper and lower meridians 
of an ogive-cylinder at M co =2.96 and a c =Q.O deg. 
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Plane of symmetry 


YRC 



Figure 17. - Convention for ordering coordinates for a 
noncircular cross section at X = XFC(J). 
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(a) Page 1 

Figure 18.- Sample case input for program DEMON3 
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Figure 18.- Continued. 
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Figure 18.- Continued. 
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Figure 18.- Continued. 
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Figure 18.- Continued. 
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Figure 18.- Continued. 
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Figure 18.- Continued. 
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Figure 20.- Sample case output from program DEMON3 
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Figure 20.- Continued. 
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inputted coefficients for numerical mapping, mnfc = 25 mxfc 
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(f) Pages 8 through 13 
Figure 20 - Continued. 
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(h) Pages 15 through 20 
Figure 20 - Continued. 
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(i) Pages 21 through 26 
Figure 20- Continued. 



body panel areas and inclination angles 

panel area delta theta delta theta 

rad rad deg deg 

1 0.23560 0.39203 -3.01069 22. 4 61 77 -172.500 
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(j) Pages 27 through 32 
Figure 20.- Continued. 



a lpha^O source strengths 

1 4 0.01303 0.01302 0.0130/] 0.01301 0.0130/1 0.01299 0.01286 0.01286 0.01299 0.01304 
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(k) Pages 32 and 33 
Figure 20.- Continued. 
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(1) Pages 33 through 38 
Figure 20.- Continued. 
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(m) Pages 38 and 39 
Figure 20.- Continued. 
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x r dr/dx m 

0000 1.1528 0.1034 1.6000 

body surface pressure distribution 
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summary of vortex field at x » 18.784 h = 1.28400 
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X r dr/dx m 

18.7840 1.8740 0.0000 1.6000 

body surface pressure distribution 
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Figure 20- Continued. 



Stratford separation criterion (laminar) f(s) =0.01511 

ummary of pressure distribution and separation points on body ... x = 18.78 

+y side: y z beta arc cp cp' dcp' /dx 

stagnation pt . 0.000 -1.952 0.000 0.019 0.020 

min. pressure 1.323 -1.114 49.907 1.652 -0.066 0.000 0.000 
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**************************************** contribution of nose section to total loads 
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(v) Pages 94 and 95 
Figure 20.- Continued. 



step 2 loads on forward-finned section with effects of forebody vorticity 
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fin section input geometry description 
fin property fin 1 or r fin 2 or 1 fin 3 or u fin 
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triplet-induced velocities at fin control points 
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vortex-induced velocities at fin control points 
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velocities from constant-u panels 
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coefficients at points on body interference shell with all perturbation velocities ** bdyshl 
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integration of the bernoulli pressure distribution ** formom ** 
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(ad) Pages 104 through 106 
Figure 20.- Continued. 



bernoulli loading 
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velocities at fin panel centroid points ** finvel ** 
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velocities at fin panel outboard aft corners ** finvel ** 
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fin loading information *** step 
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Figure 20- Continued. 



spanwise distributions 
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augmentation of fin normal force from suction conversion in proportion with factor 

kvse = 0.500 
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y/ (b/ 2 ) cn*c/{ 2 *b) ct*c/( 2 *b) cyl*c/{ 2 *b) cytot *c/ ( 2 *b) cs*c/( 2 *b) csint ybar gamnet(i) gamma, le/vinf 
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fin panel forces 
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fin vorticity due to attached flow 
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SAMPLE DEMON 3 RUN FOR NONCIRGULAR NASA MODEL: TR 2434 — DIMENSIONS IN INCHES 
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summary of vortex field at x = 33.130 h = 1.11400 
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(au) Page 130 
Figure 20 - Continued. 



Stratford separation criterion (laminar) f (s) =0.01511 
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************************************* contribution of afterbody section to total loads 
strengths and positions of vortices at end of nose section 
a - 1.5000 
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s ^ e P 4 loads on tall-finned section with effects of upstream vorticity 
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triplet -induced velocities at fin control points 
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vortex- induced velocities at fin control points 
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velocities from constant -u panels 
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coefficients at points on body interference shell with all perturbation velocities ** bdyshl 
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integration of the bernoulli pressure distribution ** formom ** 
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bernoulli loading 

rolled body-axis coordinates unrolled body-axis coordinates 
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velocities at fin panel centroid points ** finvel ** 
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velocities at fin panel outboard aft corners ** finvel 
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pressure loadings at control points 

x(j) Y(j) z(j) deltp, lin . deltp, bern. 

1.509465 1.841271 0.000000 0.169244 0.170760 
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Figure 20.- Continued. 
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